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I.  INTRODUCTION 


A.  PURPOSE  AND  OVERVIEW 

The  main  purpose  of  this  thesis  is  to  conduct  an  experiment,  using  two  different 
types  of  combat  attrition  models,  to  determine  whether  their  results  are  similar.  There 
is  a  belief  that  deterministic  Lanchester  Equation  type  models  closely  parallel  the  mean 
results  of  stochastic  computer  simulations.  This  research  will  attempt  to  determine 
whether  the  results  of  a  deterministic  "Generalized  Lanchester  Equation"  type  model 
and  a  stochastic  computer  simulation  model  are  similar.  This  work  is  the  first  such 
research  of  this  type  published  to  the  author's  knowledge. 

Chapter  One  provides  the  reader  with  an  introduction  to  modeling  and  an 
overview  of  the  history  of  Lanchester  Equation  type  models.  The  purpose  of  Chapter 
Two  is  to  introduce  the  two  models  which  are  utilized  in  this  experiment  and  Chapter 
Three  is  a  presentation  of  the  experimental  design.  The  results  of  the  experiment  are 
provided  in  Chapter  Four  and  a  discussion  of  the  conclusions  and  recommended 
follow-up  research  are  given  in  Chapter  Five. 

B.  MODELING 

Throughout  history,  military  analysts  have  been  developing  and  utilizing  various 
types  of  models  for  detailed  analysis  to  assist  in  decision  making.  A  model  is  defined 
as  a  "simplified  representation  of  the  entity  it  imitates  or  simulates"  [Ref.  1:  p.  1].  A 
second  definition  of  a  model  used  by  the  U.S.  Army  Models  Review  Committee,  "an 
abstract  representation  of  reality  which  is  used  for  the  purpose  of  prediction  and  to 
develop  understanding  about  the  real-world  processes"  [Ref.  2:  p.  5],  implies  that 
models  are  designed  to  be  close  representations  of  a  real-ivorld  entity  or  process.  A 
final  definition  for  a  military  model  is,  "an  abstraction  of  reality,  the  elements  of  which 
are  chosen  for  (a)  an  investigative  purpose  or  (b)  a  resource  management  purpose;  in 
other  words,  an  abstraction  to  assist  in  making  decisions."  [Ref.  1:  p.  3).  From  these 
definitions,  it  can  be  seen  how  important  the  development  and  use  of  models  is  to  the 
military  analyst. 

In  general,  there  are  three  types  of  models: 

1)  Iconic  model 

2)  Analogue  model 

3)  Symbolic  model 


An  iconic  model  is  a  miniature  version  of  the  entity,  such  as  an  airplane  or  a  tank 
model.  An  analogue  model  is  an  artificial  representation  of  reality.  An  example  of  an 
analogue  model  is  a  map  which  represents  the  three  dimensional  real-world  on  a  two 
dimensional,  small  scale  sheet  of  paper.  A  symbolic  model  is  one  in  which  words  or 
numerical  descriptions  are  used  to  represent  an  entity  or  process.  An  example  of  a 
symbolic  model  is  a  mathematical  equation  which  represents  a  process  such  as 
attrition. 

The  types  of  models  used  by  military  analysts  cover  the  entire  spectrum  from 
military  field  exercises  (which  can  be  thought  of  as  iconic  models)  to  analytical  models. 
Figure  1.1  is  a  combination  of  information  found  in  Taylor's  book  Force-on-Force 
Attrition  Modelling  [Ref.  2:  p.  7],  and  Hughes'  book  Military  Modeling  [Ref.  1:  p.  10], 
which  shows  the  various  types  of  models  in  use  and  their  characteristics  of  operational 
realism,  degree  of  abstraction,  convenience  and  accessibility. 


Figure  1.1  Current  Model  Types. 


The  three  most  convenient  and  accessible  types  of  models  shown  in  Figure  1.1 
are  war  games,  computer  simulations  and  analytical  models.  These  three  types  of 
models  are  also  the  most  commonly  used  to  model  combat  attrition.  A  brief 
description  of  each  of  them  will  be  given  to  better  understand  and  help  distinguish 
between  them. 

1.  War  Games 

War  games  are  conducted  using  real-world  combat  scenarios  and  people 
playing  various  positions  on  a  headquarters  battle  staff.  A  war  game  is  conducted  to 
allow  a  commander  and  his  staff  to  exercise  the  staff  decision  making  processes  in 
realistic'  combat  scenarios.  This  allows  the  commander  to  exercise  and  develop  his 
staff  without  incurring  the  true  cost  of  combat,  the  loss  of  human  life.  War  games  can 
be  run  utilizing  the  sand  table  method  or  on  computers.  The  key  to  war  gaming  is  that 
there  is  always  a  person  making  the  decisions  that  are  required  in  the  game. 

2.  Computer  Simulations 

Computer  simulations  are  also  used  to  model  the  combat  attrition  process. 
The  primary  difference  between  computer  simulations  and  war  games  is  the  method 
used  to  represent  the  decision  process.  War  games  use  people  playing  staff  roles  to 
make  the  decisions  where  simulations  use  algorithms  to  represent  the  decision  process. 
A  combat  simulation  not  only  represents  the  combat  process,  it  acts  it  out  from  start 
to  finish.  A  computer  simulation  begins  with  a  set  of  input  parameters  and  runs 
continuously  until  completion  of  the  battle  and  then  provides  the  results  of  the  battle 
as  output.  Computer  simulations  which  utilize  pseudo-random  number  generators  to 
determine  the  results  of  random  events,  such  as  the  outcome  of  one  soldier  firing  at 
another,  are  called  Monte  Carlo  Simulations  and  are  stochastic  models. 

3.  Analytical  Models 

Analytical  models  are  the  third  type  of  model  used  in  modeling  combat 
attrition.  Analytical  models  are  symbolic  models  which  use  mathematical  symbols  and 
equations  to  represent  the  combat  attrition  process.  Analytical  models  can  be 
developed  as  stochastic  models  or  deterministic  models.  Stochastic  models,  as 
introduced  in  the  previous  section,  utilize  probablity  distributions  to  determine  certain 
variables.  Therefore  the  output  from  two  consecutive  runs  of  a  stochastic  model  will 
more  than  likely  be  different.  The  deterministic  model  however,  will  produce  the  same 
set  of  output  values  for  a  given  set  of  input  parameters. 


My  research  is  focusing  on  the  comparison  of  a  stochastic  Monte  Carlo 
Simulation  and  a  deterministic  analytical  model.  The  deterministic  model  l  am  using 
fcr  comparison  is  a  Generalized  Lanchester  Equation  model  which  was  developed  from 
the  original  Lanchester  Equation  Theory.  To  provide  a  basis  for  the  readers 
understanding  of  this  experiment,  1  will  give  a  brief  introduction  to  the  history  of  the 
Lanchester  Equation  type  models. 

C.  LANCHESTER  EQUATION  THEORY 

1.  History  Of  Lanchester  Equation  Type  Modeling 

F.  W.  Lanchester  was  a  British  aeronautical  engineer,  who  in  1914  developed 
two  sets  of  simple  differential  equations  to  model  the  combat  attrition  process  between 
two  opposing  homogeneous  forces  (i.e.,  fighter  aircraft  vs.,  fighter  aircraft).  His 
intentions  were  to  provide  insight  into  the  dynamics  of  combat  under  modern 
conditions'  of  warfare  and  to  justify  the  principle  of  concentration  of  forces.  The  term 
modem  conditions  of  warfare  will  be  addressed  again  later. 

Lanchester  s  original  work  was  designed  to  model  a  force-on-force  attrition 
process  involving  two  homogeneous  forces.  He  developed  two  combat  attrition 
models,  each  having  it  s  own  unique  assumptions,  for  this  purpose.  These  original 
models  were  designed  to  model  an  aimed  lire  combat  scenario  and  an  area  lire 
combat  scenario.  These  two  models  are  often  referred  to  as  the  classical  Lanchester 
Equation  Theory. 

2.  Modeling  Aimed  Fire 

The  phrase  aimed  fire,  as  Lanchester  used  it  in  his  original  work,  refers  to 
combat  between  two  forces,  X  and  Y.  where  each  combatant  from  X  force  acquires  a 
Y  force  target  (i.e.,  locates  and  takes  aim)  and  fires.  An  example  of  this  type  combat  is 
an  infantry  battle  for  control  of  an  area  such  as  hill  224.  The  following  assumptions 
apply  for  this  model: 

1)  The  combat  being  modeled  involves  two  homogeneous  forces  (i.e..  Infantry  vs. 
Infantry). 

2)  The  entire  X  force  and  Y  force  are  within  weapons  range  of  one  another 

3)  The  effects  of  weapons  rounds  are  independent. 

4)  Each  of  the  forces  is  well  enough  aware  of  the  location  and  condition  of  all 
enemy  forces  so  that  they  will  engage  only  live  enemy  units.  Also  the  rate  at 
which  they  kill  enemy  targets  is  constant. 

5)  Fire  is  uniformly  distributed  over  surviving  enemy  targets. 


Lanchester  s  onginal  model  for  aimed  fire  is  given  by  Equation  Set  1.1: 


dx  dt  -  -ay(t) 

with  x(0) 

dy  dt  -  -bx(t) 

with  y(0) 

(eqn  1.1 ) 


where  x(t)  is  defined  as  the  number  of  combatants  of  X  force  at  time  t,  y(t)  is  defined 
as  the  number  of  combatants  of  Y  force  at  time  t.  and  a  and  b  are  constants  which  are 
called  attntion-rate  coefficients.  The  attrition-rate  coefficient  a  represents  the 
effectiveness  of  a  Y  combatant  killing  a  X  combatant  per  unit  of  time.  Similarly  the 
attntion-rate  coefficient  b'  represents  the  effectiveness  of  a  X  combatant  killing  a  Y 
combatant  per  unit  of  time.  This  set  of  equations  assumes  that  the  rate  at  which  a 
force  X  can  cause  casualties  to  a  Y  force  is  proportional  to  the  number  of  combatants 
in  the  X  force.  Similarly,  the  rate  at  which  a  combatant  from  Y  force  can  cause 
casualties  to  the  X  force  is  assumed  to  be  proportional  to  the  number  of  combatants  in 
the  Y  force.  When  Equation  Set  1.1  is  integrated,  the  resulting  state  equation  is 
obtained  and  has  been  labeled  the  Lanchester  Square  Law; 

b(x02-x(t)2)«a(y02-y(t)2).  (eqn  1.2) 

3.  Modeling  Area  Fire 

The  phrase  area  fire  as  Lanchester  used  it  in  his  onginal  work,  refers  to  a 
combat  scenano  involving  two  homogeneous  forces,  say  X  and  Y,  which  are  uncertain 
of  exact  enemy  locations  and  engage  one  another  by  firing  in  the  general  area  where 
the  enemy  force  is  located.  This  scenano  has  the  following  assumption  set: 

1)  The  combat  being  modeled  involves  two  homogeneous  forces  (i.e..  Infantry  vs. 
Infantry) 

2)  The  entire  X  force  and  Y  force  are  within  weapons  range  of  one  another. 

3)  The  effects  of  weapons  rounds  are  independent. 

4)  Each  force  is  aware  only  of  the  general  location  of  the  enemy  force  and 
therefore  engages  the  enemy  by  firing  into  that  general  area  without  the 
benefit  of  knowing  their  effectiveness. 

5)  Fire  from  all  surviving  combatants  is  uniformly  distributed  over  the  area 
which  the  enemy  occupies. 

6)  Each  force  has  the  same  vulnerable  area  to  enemy  lire. 


Lanchesters  original  model  for  area  fire  is  defined  by  the  Equation  Set  1.3: 


dx  dt  -  -axy, 
dy'dt  ”  -byx. 


(eqn  1.3) 


This  model  assumes  that  the  rate  of  attrition  is  proportional  to  both  force  levels  x(t) 
and  y(t).  This  means  that  in  the  area  fire  model,  the  rate  at  which  a  X  force  can  cause 
casualties  to  a  Y  force  is  not  only  dependant  upon  the  number  of  combatants  in  its 
own  force,  but  also  on  the  number  of  combatants  alive  in  the  Y  force.  The  rate  at 
which  the  Y  force  can  cause  casualties  on  the  X  force  is  similarly  dependent  on  both 
their  own  force  level  and  the  force  level  of  the  opposing  force  X. 

The  well  known  'Linear  Law'  is  obtained  by  integration  of  Equation  Set  1.3 
and  is  shown  here  in  Equation  1.4  : 

b(x0-x(t))«a(y0-y(t»),  (eqn  1.4) 

where  xQ  and  yQ  are  the  initial  force  levels  at  time  t-0.  x(t)  and  y(t)  are  the  force 
levels  at  time  t,  and  a  and  b  are  the  attrition-rate  coefficients  for  each  force. 

4.  Attrition  Rate  Coefficients 

The  Lanchester  attrition-rate  coefficients  were  introduced  in  the  models  above 
and  were  defined  as  constants.  In  order  to  use  the  Lanchester  models  there  must  exist 
a  set  of  acceptable  attrition-rate  coefficients  for  each  possible  type  of  engagement. 
Lanchester  s  original  model  utilizing  homogeneous  forces,  only  requires  the  availability 
of  a  limited  number  of  coefficients.  However,  in  Lanchester  Equation  type  models 
which  allow  for  the  combat  between  heterogeneous  forces,  the  number  of  attrition-rate 
coefficients  grows  very  rapidly  when  one  considers  all  the  possible  combinations  of 
combatant  weapons  vs.  combatant  weapons  that  can  occur.  As  Lanchester  postulated 
his  original  work,  he  used  constants  for  attrition-rate  coefficients,  which  implies  that 
the  kill  capability  of  a  force  does  not  change  over  time.  With  constant  attrition-rate 
coefficients  these  models  are  easy  to  solve.  However,  there  is  no  reason  why  we  should 
believe  that  such  factors  as  range,  weather,  visibility,  and  training  do  not  affect  the  kill 
capability  of  a  combatant.  Taylor  provides  two  methods  for  determining  numerical 
values  for  Lanchester  attntion-rate  coefficients  which  are  used  in  the  United  States. 


1)  A  statistical  estimate  based  on  combat  data  generated  by  a  detailed  Monte 

Carlo  combat  simulation. 

2)  An  analytical  submodel  of  the  attrition  process  for  the  particular  combination 

of  firer  and  targets. 

In  the  first  method,  the  simulation  output  is  used  to  fit  one  or  more  free 
parameters  in  the  analytical  model  in  an  attempt  to  have  it  to  provide  similar  results. 
S.  Bonder  has  labeled  this  approach  a  'fitted-parameter  analytical  model  [Refs.  2,3:  pp. 
45.73-88). 

The  second  method  for  developing  the  attrition-rate  coefficients  has  been 
labeled  an  'independant  analytical  model'  by  S.  Bonder  [Refs.  2,3:  pp.  47,73-88],  The 
concept  for  developing  the  coefficients  is  to  consider  a  single  firer  shooting  at  a  single 
target  (which  is  stationary  and  does  not  fire  back).  These  coefficients  are  developed 
under  perfect  conditions  with  this  method.  There  has  been  a  large  amount  of  work 
done  in  the  development  of  acceptable  attrition-rate  coefficients.  Taylor  provides  a 
good  list  of  references  in  his  book  Force-on-Force  Attrition  Modeling  [Ref.  2). 

5.  Modern  Conditions  of  Warfare 

In  1914  when  F.  W.  Lanchester  developed  his  original  models,  he  was 
attempting  to  model  the  combat  attrition  process  so  that  he  could  learn  more  about 
the  dynamics  of  combat  under  modem  conditions  .  The  phrase  modem  conditions 
was  an  important  one  then  as  it  is  today.  In  earlier  sections  Lanchesters  models  were 
introduced  and  the  assumptions  that  applied  to  them  were  given.  The  assumption  sets 
are  what  tailor  the  model  to  the  conditions.  Some  of  these  original  assumptions  no 
longer  seem  appropriate  to  model  what  we  would  term  combat  under  modern 
conditions'.  The  next  section  discusses  some  of  the  shortcomings  of  the  original 
Lanchester  Models,  but  before  that  it  is  important  to  lay  out  a  solid  definition  for 
combat  under  modem  conditions. 

Today  the  conditions  on  the  combat  battlefield  are  very  complex.  First,  the 
air  and  land  battles  are  thought  of  as  one  battle  and  given  the  title  AIR-LAND 
BATTLE'.  The  ground  forces  commander  from  the  battalion  level  on  up.  has  a 
combination  of  different  combat  units  under  his  direct  command.  An  Army  Division 
Commander  has  Infantry,  Armor,  and  Artillery  Battalions  all  assigned  to  his  command 
The  battalions  which  are  assigned  are  not  pure  units,  they  too  have  a  combination  of 
combat  units  assigned.  This  concept  has  been  adopted  by  the  Army  as  the  way  to 
fight  the  next  war.  The  Army  has  orgam/cd  under  the  combined  arms  concept  and 
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believe  that  placing  a  mix  of  supporting  combat  forces  together  on  the  battlefield  will 
be  a  force  multiplier. 

The  Navy  has  developed  much  the  same  type  of  conditions  at  sea.  The 
Carrier  Battle  Groups  (CBG)  are  combinations  of  many  types  of  naval  combat  power 
and  therefore  brings  a  variety  of  forces  together  to  fight  a  battle  at  one  point  under  the 
Ofiicer-in-Tactical-Command  (OTC). 

6.  Shortcomings  With  The  Original  Lanchester  Models 

Lanchesters  original  models  provided  a  very  simplified  model  of  modem 
warfare  which  has  served  as  the  basis  for  a  great  deal  of  further  work  in  this  area. 
There  are  many  shortcomings  of  the  original  models  which  have  been  identified  and 
have  stimulated  much  of  the  work  done  in  the  years  following  their  development. 
Lanchesters  models  have  been  popular  because  of  their  simplicity  and  the  ease  of 
tracing  through  the  mathematical  computations.  In  this  section  I  will  discuss  some  of 
the  shortcoming  which  have  been  identified  in  applying  Lanchesters  original  laws  to 
model  today  s  combat. 

First,  not  all  of  today  s  combat  scenarios  can  be  neatly  placed  into  the 
square'  or  linear  attrition  models.  The  most  often  referenced  and  easily  explained 
example  is  a  guerilla  ambush  scenario.  In  this  case  it  is  assumed  that  the  guerilla  force 
is  heavily  camouflaged  and  the  force  being  attacked  is  in  the  open.  It  is  apparent  that 
the  attacked  force  would  suffer  a  squared'  attrition  rate  from  the  aimed  fire  of  the 
guerilla  force  which  has  clear  line  of  sight  of  the  entire  force  under  attack.  On  the 
other  hand,  the  defending  force  would  not  have  clear  line  of  sight  of  the  attacking 
guerilla  force  and  therefore  would  not  be  properly  modeled  by  a  square  attrition 
model,  but  rather  the  linear  attrition  model.  However,  even  this  does  not  properly 
model  the  scenario.  With  time  the  defending  force  will  be  able  to  locate  the  attacking 
force  well  enough  that  their  fire  should  be  considered  as  aimed  fire  and  modeled  by  a 
squared  attrition  model.  This  illustrates  the  dynamic  conditions  on  the  battlefield 
which  must  be  considered.  [Ref.  4] 

One  of  the  main  assumptions  of  Lanchester  s  laws  is  that  the  two  opposing 
forces  are  comprised  of  homogeneous  units  lie.  Infantry  vs.  Infantry)  This  is  rarely 
the  case  on  the  battlefield  today  Today  the  doctrine  for  a  land  battle  is  the  Air-Land 
Battle  Doctrine,  and  is  written  for  combined  arms  units  1  he  combined  arms  concept 
was  introduced  earlier,  but  warrants  reemphasis  I  lie  force  that  must  be  modeled 
today  combines  several  different  combat  units  under  one  commander  to  light  a  battle 


An  additional  consideration  is  the  close  relationship  between  the  ground  force  and  the 
close  air  support  provided  by  the  Air  Force  which  adds  another  unit  type  into  the 
model.  The  same  type  of  unit  composition  is  currently  found  in  the  Navy  Carrier 
Battle  Groups.  The  need  for  a  model  which  allows  for  several  different  types  of  units 
(the  heterogeneous  force  model)  is  obvious.  The  use  of  heterogeneous  forces  requires 
the  model  to  consider  allocating  fire  power  on  a  variety  of  targets.  The  factors  that 
enter  in  the  allocation  process  are;  the  attrition  rate  of  each  type  of  weapon  for  each 
of  the  opposing  forces  weapons  types,  the  number  and  type  of  enemy  targets,  and  the 
relative  significance  of  each  of  the  enemy  targets. 

Another  assumption  of  the  original  Lanchester  theory  is  that  all  of  the  forces 
of  both  sides  are  committed  to  battle  at  the  beginning  and  there  are  no  reinforcements 
available.  Furthermore,  there  is  no  allowance  for  the  possibility  of  withdrawing  some 
or  all  of  a  unit's  forces.  These  factors  are  definitly  consider  by  the  commander  on  the 
battlefield  today. 

Separation  of  forces  on  the  battlefield  is  ignored  in  the  original  models.  They 
assume  that  the  units  remain  fixed  in  location  and  are  always  within  range  of  the 
opposing  forces  weapons.  To  better  model  the  combat  process,  a  movement  factor 
which  would  affect  the  weapons'  ranges  and  the  accuracy  of  each  weapon  on  the 
battlefield  would  be  required. 

The  original  model  is  designed  to  model  a  'fight  to  the  finish'  type  of  combat. 
This  definition  for  battle  termination,  of  fighting  to  the  last  man,  is  not  realistic  on 
today  s  battlefield.  There  is  some  point  in  a  battle  where  the  defender  will  decide  to 
withdraw  or  surrender  before  the  entire  force  is  destroyed.  [Ref.  4] 

With  this  introduction  to  Lanchester  Equation  type  models  as  a  basis,  the 
next  chapter  will  present  the  two  models  used  for  this  research. 


II.  THE  MODELS 


A.  GENERAL 

Two  models  were  used  for  this  research,  one  is  a  deterministic  model  and  the 
other  is  a  stochastic  computer  simulation  model.  The  deterministic  model  used  was 
developed  at  the  Naval  Postgraduate  School  bv  Professor  Paul  H.  Moose  and 
Professor  John  M.  Wozencraft.  Their  model  is  a  Generalized  Lanchester  Equation 
type  model.  It  will  be  referenced  throughout  the  remainder  of  this  thesis  as  the  M  W 
model  or  the  Moose  Wozencraft  model. 

The  stochastic  computer  simulation  model  which  was  utilized  for  this  research  is 
a  Monte  Carlo  Simulation  written  by  Professor  Don  E.  Harrison,  Jr.  at  the  Naval 
Postgraduate  School.  It  will  be  referenced  throughout  the  remainder  of  the  thesis  as 
the  C  S  model  or  the  computer  simulation  model. 

In  the  following  sections,  each  of  these  models  will  be  presented  to  provide  a 
general  understanding  of  there  characteristics. 

B.  MOOSE/WOZENCRAFT  MODEL 

Paul  H.  Moose  is  a  professor  in  the  Electrical  Engineering  Department  at  the 
Naval  Postgraduate  School  and  John  M.  Wozencraft  was  a  professor  with  the 
Electrical  Engineering  Department  prior  to  his  retirement.  Both  have  been  associated 
with  the  academic  group  for  the  Joint  Command,  Control,  and  Communications 
Curriculum,  and  have  been  interested  in  the  decision  and  control  problems  a  modern 
day  military  commander  faces.  They  wanted  to  develop  a  better  understanding  of  how 
command  and  control  decisions  affect  battle  outcome.  To  do  this,  an  analytical  model 
of  the  attrition  process  that  would  provide  an  adequate  representation  of  the  modem 
conditions  of  combat  was  needed.  Such  a  model  required  modeling  a  variety  of  forces 
as  an  aggregate  force  fighting  on  the  battlefield.  [Ref.  5] 

As  presented  in  chapter  one,  the  military  commander  in  today's  environment  has 
several  types  of  units  available  to  him  to  fight  a  battle.  The  question  of  optimum 
resource  allocation  involves  the  use  of  a  variety  of  these  units.  For  the  model  to  be 
useful  it  had  to  be  easy  to  interpret  and  understand.  Lanchester  Equation  type  models 
were  easy  to  interpret,  but  were  limited  in  their  ability  to  model  a  variety  of  forces. 
This  motivated  their  research  into  the  dynamics  of  a  Generalized  System  of  Lanchester 
Equations  of  the  type  shown  here  in  Equation  Set  2.1: 


(eqn  2.1) 


Xj(t)=  F^Xj.y.t): 


yj(t)=  Gj(x,yj,t): 


where  i=  1,2 . N, 


where  j=  1,2 . M, 


where  N  components  of  a  non-homogeneous  X-force  engage  M  components  of  a  Y- 
force.  As  shown  in  Equation  Set  2.1,  the  number  of  combatants  of  type  xj  is  a 
function  of  their  own  forces,  the  total  Y  force  (the  sum  of  the  yj  forces)  available  to 
fight,  and  time.  Similarly,  the  rate  of  change  of  a  yj  force  is  a  function  of  it  s  own 
force  level,  the  total  X  force  (the  sum  of  all  xj  forces)  available  and  time. 

The  following  set  of  equations  represent  the  Moose/ Wozencraft  model: 


xi(t)  =  -uixi(t)-JaijXi(t)yj(t)-JbijVj(t)  +  rj  where  i=  1,2 . X, 

J  -J 

yj(t)“-Vjyj(t)-JcijXi(t)yj(t)-XdijXi(t)  4-  Sj  where  j  *  1,2 . M, 


(eqn  2.2) 


where 

1)  The  Xj  and  yj  represent  Blue  and  Orange  forces  of  different  types 

2)  The  u-  and  Vj  are  self  attrition  coefficients 

3)  The  a-  and  c-  are  area  fire  coefficients 

4)  The  b-j  and  d-  are  aimed  fire  coefficients 

5)  The  r-  and  Sj  are  resupply  coefficients 

6)  Xj.Vj  £  0  for  all  i  and  j 

7)  £xrX,£  yj-Y. 

*  J 

The  M/W  model  is  designed  as  an  NxM,  heterogeneous  force,  combat  attrition 
model.  The  significance  of  this  fact  is  that  it  allows  for  closer  modeling  of  the  modem 
combat  force  discussed  in  Chapter  One. 

In  the  M/W  model,  Equation  Set  2.2,  several  things  are  different  from  the 
original  Lanchester  Equation  models  presented  in  Chapter  One.  One  major  difference 
is  that  the  M/W  model  allows  for  modeling  N  type  X  forces  and  M  type  Y  forces.  The 
Moose/Wozencraft  model  treats  the  attrition  of  a  given  force  as  one  process  where 
Lanchester's  original  work  broke  the  attrition  process  into  two  models,  aimed  and  area 
fire  models.  With  that  fact  in  mind  then,  the  rate  of  change  of  a  force  x-  is  a  function 
of  four  separate  factors;  self-attrition,  area  fire  attrition,  aimed  fire  attrition  and 
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resupply.  In  the  case  of  determining  the  attrition  to  force  Xj  ,  the  self-attrition  is 
modeled  using  a  self-attrition  coefficient  -Uj  which  is  multiplied  by  the  current  Xj  force 
level.  This  yields  the  attrition  to  force  xj  w’hich  is  attributed  to  'self  losses'  such  as 
disease  and  defections.  The  area  fire  attrition  is  modeled  similarly  to  Lanchesters 
original  model,  with  the  added  capability  to  account  for  all  the  possible  combinations 
of  y:  forces  which  can  cause  casualties  to  x ■  by  means  of  area  fire.  Thus  an  attrition- 
rate  coefficient  for  each  possible  combination  of  forces  engaging  in  battle  must  be 
available.  Each  Vj  force  has  an  area  fire  coefficients  -ay  associated  with  it.  The  sum  of 
the  attrition-rate  coefficients  (a-j)  multiplied  by  the  Vj  force  level  can  then  be  thought 
of  as  an  aggregate  area  fire  attrition  rate  which  multiplies  the  x-  force  level.  The  aimed 
fire  is  very  similar.  It  is  also  similar  to  the  original  Lanchester  model  with  the  added 
capabilities  to  model  NxM  force  battles.  The  aimed  fire  is  modeled  by  multipling  each 
>'j  force  by  an  attrition-rate  coefficient  -by  and  summing  them  up.  This  is  then  the 
aggregate  aimed  fire  attrition  rate.  Finally  the  last  factor  modeled  is  the  resupply  of 
each  Xj  force  with  a  resupply  variable  rj.  The  sum  of  the  self-losses,  area  fire  losses, 
aimed  fire  losses  and  resupply  is  equal  to  the  change  in  the  x-  force. 

The  model  determines  the  change  in  a  Vj  force  in  exactly  the  same  manner.  The 
attrition  rate  coefficients  are  different  of  course,  but  the  process  is  the  same.  The  self¬ 
attrition  coefficient  for  Vj  is  -Vj.  The  area-fire  attrition  coefficients  are  given  as  -Cy  and 
the  aimed  fire  coefficients  are  given  as  -dy.  The  resupply  variable  for  yj  is  Sj. 

The  relative  usefulness  of  this  model  for  modeling  the  combat  attrition  process 
on  today's  battlefield  is  apparent  from  this  discussion.  First,  it  provides  the  capability 
to  model  a  battle  involving  two  heterogeneous  forces  which  we  saw  is  necessary  to 
analyze  the  optimum  force  allocation  issue.  Secondly,  it  is  a  deterministic  model,  so  it 
is  easy  to  understand  and  interpret.  It  has  many  similarities  to  the  original  Lanchester 
Equation  models  which  are  widely  studied  and  understood. 

The  results  of  Professor  Moose  and  Professor  Wozencraft's  research  into  the 
dynamical  properties  of  this  model  are  provided  in  a  paper  which  they  are  submitting 
to  the  Military  Operations  Research  Society  for  publication  titled  Characteristic 
Trajectories  of  Generalized  Lanchester  Equations  [Ref.  6]. 

C.  STOCHASTIC  COMPUTER  SIMULATION  MODEL 

The  stochastic  model  is  a  Monte  Carlo  Simulation  designed  to  test  the 
Moose/ Wozencraft  model.  The  basis  of  this  model  was  designed  by  Don  E.  Harrison, 


Jr.,  a  professor  at  the  \aval  Postgraduate  School.  Professor  Harrison  was  asked  to 
design  a  stochastic  model  of  the  attrition  process  which  included  the  same  factors  as 
the  Moose, 'Wozencraft  model.  The  combination  of  his  efforts  and  some  additional 
code  to  run  the  program  for  a  given  number  of  replications  and  calculate  statistical 
data  resulted  in  the  model  shown  in  Appendix  B  entitled  Computer  Simulation  model 
(C/S). 

The  C/S  model  will  be  presented  in  the  following  sections: 

1)  Input  parameters 

2)  The  Combat  Cycle 

3)  The  Output 

The  material  in  each  section  is  presented  in  a  condensed  users  manual  form  to  provide 
the  reader  with  an  understanding  of  the  model  and  the  capability  to  exercise  the  model 
if  desired. 

I.  Input  parameters 

The  computer  code  is  written  in  FORTRAN77  and  uses  an  exec  file  to  define 
all  the  input  files,  output  files,  load  and  execute  the  program. 

The  input  file  contains  all  the  required  input  data.  Figure  2.1  shows  the 
contents  and  format  of  the  input  data  file.  The  first  line  of  the  file  is  a  heading  which 
is  read  by  the  program  and  allows  the  user  to  label  and  keep  track  of  which  test  data 
set  is  being  utilized.  Next  the  data  required  for  the  Xj  forces  is  given.  On  the  second 
line  the  number  of  types  of  X  forces  is  given  (the  va  iable  NTYPX  in  the  program), 
this  corresponds  to  the  N  value  in  the  M'W  model  which  specifies  the  number  of  types 
of  X  forces.  Then  the  initial  force  level  for  each  Xj  force  is  given,  their  locations,  their 
corresponding  self-attrition  coefficients  (u-)  and  the  resupply  variables  (r-).  Next,  the 
attrition-rate  coefficients  are  provided  in  matrix  form.  The  aimed  fire  attrition-rate 
coefficients  'd-'  are  given  in  lines  three  through  five.  Area  fire  attrition- rate 
coefficients  'cjj'  are  given  in  lines  six  through  eight.  Then  the  required  information  for 
the  yj  forces  is  listed.  Once  again,  the  number  of  type  of  yj  forces  is  given  (the 
N’TYPY  variable  in  the  program)  which  corresponds  to  M  in  the  M/W  model.  Then 
each  Vj  initial  force  level,  location,  self  attrition-rate  coefficients  Vj'  and  resupply 
variables  'Sj'  are  given.  Lines  ten  through  twelve  contain  the  aimed  fire  coefficients 
bjj'  and  lines  thirteen  through  fifteen  contain  the  area  fire  coefficients  'a-'.  Line 
sixteen  of  the  file  contains  seed  numbers  for  the  pseudo-random  number  generators. 


1  01/16/87  TEST  OF  (2X2),  1  UNSTABLE  ROOT  S 

2  NSTEP 

3  NTYPX  NTX1  NTX2  NTX3  loci  loc2  loc3  Uj  u2  u3  r,  r2  r3 

4  dj j  dj 2  d,3 

5  d21  d22  d23 

6  d31  d32  d33 

7  C11  c12  c13 

8  C21  c22  c23 

9  C31  c32  c33 

10  NTYPY  N'TYl  NTY2  NTY3  loci  loc2  loc3  v2  v3  Sj  s2  s3 

11  bll  b12  b13 

12  b21  b22  b23 

13  b31  b32  b33 

14  all  a12  a13 

15  a21  a22  a23 

16  a31  a32  a33 

17  DSEED  DSEED  DSEED 


Figure  2.1  Input  Data  File. 

The  program  is  designed  to  be  interactive  before  it  begins  the  combat  cycle 
which  allows  the  user  to  change  the  heading  of  the  input  file  if  so  desired.  This  also 
enables  the  user  to  ensure  that  the  input  file  is  the  one  desired. 

2.  The  Combat  Cycle 

The  combat  cycle  is  a  stochastic  process.  During  the  combat  cycle  the 
combatants  are  chosen  at  random  and  the  probability  of  self-loss,  aimed  fire  loss,  and 
area  fire  losses  are  all  tested  by  use  of  a  pseudo-random  number  between  0.0  and  1.0. 
The  cycle  begins  by  randomly  selecting  which  force,  X  or  Y,  will  fire.  Once  the  force 
has  been  selected  then  a  combatant  from  the  force  is  selected.  For  the  purpose  of 
presentation,  say  that  Xj(l)  (that  is  combatant  one  of  force  type  Xj)  is  selected  to  be 
the  firer,  then  the  first  thing  that  occurs  in  the  cycle  is  a  check  to  see  if  combatant 
Xj(l)  is  a  self  loss.  This  is  accomplished  by  use  of  a  random  number  (between  0.0  and 


1.0)  which  is  compared  against  the  appropriate  self  attrition  coefficient  Uj.  If  the 
combatant  is  not  a  self  loss  then  a  Y  force  combatant  is  selected  at  random  as  a  target, 
say  >’2(2)  (the  second  combatant  of  force  type  y2)  as  an  example.  Once  the  firer  and 
target  have  been  identified,  the  next  step  tests  the  effectiveness  of  the  shot'  by  using  a 
random  number  to  test  whether  the  target  is  killed.  This  process  involves  first 
determining  if  the  target  is  killed,  which  is  done  by  testing  whether  a  random  number  is 
less  than  the  total  kill  probablity  (which  is  the  sum  of  the  aimed  fire  and  area  fire 
attrition-rate  coefficients,  or  their  kill  probabilities).  If  the  random  number  is  less,  then 
the  target  has  been  killed  and  if  not  the  shot  missed.  If  the  target  was  killed  then  this 
same  process  is  followed  to  determine  whether  it  was  a  loss  due  to  aimed  fire  or  area 
fire.  Required  data  is  tabulated  on  the  number  of  combatants  of  type  y:  which  are 
killed  by  type  xj  forces  and  each  force  is  resupplied.  The  combat  cycle  is  repeated  in 
this  manner  until  all  combatants  have  had  an  opportunity  to  fire  during  each  timestep. 
This  process  is  continued  for  each  timestep  until  one  of  two  conditions  occurs: 

1)  One  force  reaches  a  break  point  which  has  been  defined  as  NXSTOP  and 
NYSTOP.  This  simulates  the  level  of  attrition  a  commander  will  suffer  before 
pulling  back  or  surrendering  rather  than  fighting  to  the  last  man. 

2)  The  combat  cycle  has  been  repeated  for  a  given  number  of  timesteps,  which  is 
specified  in  the  input  data  file. 

3.  The  Output 

The  C/S  model  was  designed  to  provide  several  output  files  for  data  analysis. 
The  model  generates  two  files  which  provide  a  very  detailed  timestep-by-timestep 
recording  of  the  combat  cycle  outcome.  These  files  require  a  great  deal  of  storage 
space  and  therefore  during  the  interactive  portion  of  the  program  the  user  is  asked  if 
they  want  these  files  printed  in  their  complete  form.  It  is  not  recommended  that  they 
be  printed  if  the  user  desires  multiple  replications  because  the  storage  space  required  is 
too  large.  One  of  these  files  provides  a  detailed  listing  of  each  timestep  results 
including  a  summary  of  the  force  levels,  the  number  of  losses  by  each  type  of  fire  and 
the  number  of  combatants  that  were  resupplied.  An  example  of  this  file's  output  is 
shown  in  Figure  2.2.  The  other  file  contains  similar  listings  with  the  number  of  dead 
combatants  of  each  type  provided  as  well.  An  example  of  this  file  s  output  is  shown  in 
Figure  2.3. 

The  remainder  of  the  files  are  designed  to  provide  the  force  levels  of  each  force 
(i.e.,  Xj  and  yj)  for  each  timestep.  This  provides  the  data  which  can  be  used  to  plot 


02/02/87  TEST  OF  13X3), 

9  X  TYPES  NUMBER  EACH  LOCATIONS  SELF-LOSS  RE-SUPPLY 

3  200  150  175  1.0  1.0  1.0  0.005  0.010  0.020  4  4  11 

KILL  PROBABILITIES  FOR  FORCE  X(I)  SHOOTS  YIJ) 

PK(  1*1)  PK(  1»2  )  PK(  1  >3  )  PK(  2,1 )  PK(  2,2 )  PK<2,3)  PK(3,1)  PK<3,2)  PK(3,3) 
0.06000  0.02000  0.02000  0.00000  0.00000  0.00000  0.03000  0.01000  0.02000 


AREA  FIRE  COEFFICIENTS  FORCE  X(  I )  SHOOTS  YU) 

CC(1»1)  CCll, 2)  CCll,  3)  CC(  2,1 )  CC(  2,2 )  CC(2,3)  CC(3,1)  CC(3,2)  CC(3,3) 

0.00000  0.00000  0.00000  0.00000  0.00030  0.00020  0.00000  0.00000  0.00010 

9  Y  TYPES  NUTBER  EACH  LOCATIONS  SELF-LOSS  RE-SUPPLY 

3  150  125  225  1.0  1.0  1.0  0.003  0.020  0.030  4  5  14 

KILL  PROBABILITIES  FOR  FORCE  Y(I)  SHOOTS  X(J) 

PKl  1,1)  PK 1 1 , 2  )  PKI1.3)  PKl  2,1 )  PK(2,2)  PK(2,3)  PK(3,1)  PKI3,2)  PK(3,3) 

0.05000  0.01000  0.05000  0.00000  0.00000  0.00000  0.02000  0.01000  0.03000 


AREA  FIRE  COEFFICIENTS  FORCE  Y(I)  SHOOTS  X(J) 

AAt  1,1 )  AAl  1,2  )  AA(  1,3  )  AA(2,1)  AA(2,2)  AA(2,3)  AA(3,1)  AA(3,2)  AAI3,3) 
0.00000  0.00000  0.00000  0.00000  0.00020  0.00030  0.00000  0.00000  0.00010 

DOUBLE  PRECISION  RANDOM  VARIANT  SEEDS  INTEGER  SEED  VALUES 

0 . 345215790+09  0.4563576 10 *08  0.893427610*08  1209308051  1194034571  1196766258 


9  NTXS 

0  525 

NTX  NX1  NX2 
525  200  150 

NX3  RX1  RX2  RX3  NTYS  NTY  NY1 
175  4  4  11  500  500  150 

NY2  NY3  RY1  RY2  RY3 
125  225  4  5  14 

TIMESTEP 

1 

AIMED  FIRE  KILLS  OF  Y"S  BY  X  FORCE 

0(1,1) 

0(1,2)  0(1,3) 

3  0 

0(2,1)  0(2,2)  0(2,3)  0(3,1) 

2  0  0  0 

0(3,2) 

1 

0(3,3) 

1  0 

AIMED  FIRE  KILLS  OF  X"S  BY  Y  FORCE 

B(  1,1 ) 

B(  1,2  )  B  ( 1,3) 

2  1 

B(  2,1 )  B  (  2,2)  B(  2 ,3  )  B(3,l) 

10  0  0 

B( 3,2  ) 

2 

Bl 3,3 ) 

0  2 

AREA  FIRE 

KILLS  OF  Y"S  BY  X  FORCE 

Cll.l) 

C( 1,2 )  C( 1,3 ) 
0  0 

C( 2,1)  C( 2 ,2  )  C( 2,3 )  Cl 3,1 ) 
0  0  5  7 

C( 3,2 ) 

0 

C( 3,3 ) 

0  2 

AREA  FIRE 

KILLS  OF  X"S  BY  Y  FORCE 

A(  1,1) 

A(  1,2 )  All, 3) 
0  0 

A(  2 ,1 )  A(  2 ,2 >  A(  2 ,3  )  A(3,l) 

0  0  2  5 

A(  3,2 ) 

0 

A( 3,3 ) 

0  3 

9  NTXS 

1  525 

NTX  NX1  NX2 
525  200  150 

NX3  0X1  0X2  0X3  NTYS  NTY  NY1 
175  5  5  13  500  500  150 

NY2  NY3  0Y1  0Y2  DY3 
125  225  4  9  20 

4  4  XI 


TOTAL  X  FORCE  LOSS  HAS  23 >  RESUPPLY  HAS: 

AIMED  FIRE  LOSSES:  413 
AREA  FIRE  LOSSES:  028 

SELF  LOSSES:  122 
TOTAL  Y  FORCE  LOSS  HAS  33,  RESUPPLY  HAS:  4  5  14 

AIMED  FIRE  LOSSES:  412 
AREA  FIRE  LOSSES:  059 

SELF  LOSSES:  039 

♦  1  525  521  199  149  173  4  4  11  500  490  150  121  219  4  5  14 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  Y  FORCE 

AA(  1,1 )  AA(  1,2 )  AA(  1,3 )  AA(2,1)  AA(2,2)  AA(2,3)  AAI3,1)  AAI3,2)  AAI3.S) 
0.00000  0.00000  0.00000  0.00000  0.00011  0.00023  0.00000  0.00000  0.00008 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  Y  FORCE 

BBC  1,1)  88(1,2)  88(1,3)  88(2,1)  BB(2,2)  B8(2,3)  BB(3,1>  88(3,2)  88(3,3) 
0.03500  0.02333  0.02000  0.00000  0.00000  0.00000  0.02333  0.00000  0.02667 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE 

CCI1,1)  CCll, 2)  CC(  1,3)  CC( 2,1 )  CC(2,2)  CC(2,3)  CC(3,1)  CC(3,2)  CC(3,3) 
0.00000  0.00000  0.00000  0.00000  0.00027  0.00021  0.00000  0.00000  0.00005 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE 

00(1,1)  DD(  1,2)  00(1,3)  00(2,1)  00(2,2)  DD(2,3)  0D(3,1)  00(3,2)  00(3,3) 
0.05000  0.00000  0.02222  0.00000  0.00000  0.00000  0.01905  0.02286  0.00000 

STOPPED  BY  PROGRAM  AT  NSTOP  =  1 


Figure  2.2  Example  File  6  Output,  (cont'd.) 

force  level  trajectories,  which  allows  for  ease  in  performing  a  comparison  of  the  two 
model  s  results  for  a  given  input  set. 

Now  that  each  of  the  models  has  been  presented,  chapter  three  will  discuss  the 
design  of  the  experiment,  and  provide  samples  of  the  output  files  used  for  plotting  the 
force  level  trajectories. 


02/02/87  TEST  Of  <3X31, 

»  MTXS  NTX  tea  NX2  MX3  DX1  0X2  0X3  NTYS  NTY  KYI  MY 2  NY  3  0Y1  DY2  DY J 
0  525  525  200  150  175  0  0  0  500  500  150  125  225  0  0  0 

D,B,C,A  BY  RONS  IN  OESCENOINB  OROER 

1302000110 
1211000202 
1000057002 
1000025003 
TOTAL  X  FORCE  LOSS  NAS  23,  RESUPPLY  MAS:  4  4  11 

AIMED  FIRE  LOSSES:  413 

AREA  FIRE  LOSSES:  0  2  B 

SELF  LOSSES:  122 

TOTAL  Y  FORCE  LOSS  MAS  33,  RESUPPLY  MAS:  4  5  14 

AIMED  FIRE  LOSSES:  412 
AREA  FIRE  LOSSES:  054 

SELF  LOSSES:  034 

1  525  521  144  144  173  5  5  13  500  440  150  121  214  4  4  20 

A ,B ,C ,0  BY  RONS  IN  OESCENOINB  ORDER 

AV  0.00000  0.00000  0.00000  0.00000  0.00011  0.00023  0.00000  0.00000  0 . 00000 

AV  0.03500  0.02333  0.02000  0.00000  0.00000  0.00000  0.02333  0.00000  0.02447 

AV  0.00000  0.00000  0.00000  0.00000  0.00027  0.00021  0.00000  0.00000  0.00005 

AV  0.05000  0.00000  0.02222  0.00000  0.00000  0.00000  0.01405  0.02204  0.00000 

STOPPED  BY  PROGRAM  AT  NSTOP  »  1 


III.  RESEARCH  METHODOLOGY 


A.  OBJECTIVE 

The  primary  objective  of  this  research  is  to  determine  whether  a  deterministic 
combat  attrition  model  is  a  good  approximation  for  the  mean  result  of  ^  stochastic 
model  of  the  combat  attrition  process  With  that  objective  in  nund.  an  experiment  was 
designed  to  preform  a  comparison  of  two  such  models  The  models  used  in  this 
experiment  are  the  Moose  Wozencraft  model  iM  W  model),  which  is  a  deterministic 
model,  and  a  Monte  Carlo  Simulation  (<  S  model)  which  were  introduced  m  chapter 
two  A  secondary  objective  is  to  see  if  the  two  models  behave  with  similar 
characteristics  as  the  results  of  Professor  Moose  and  Professor  Wozencraft  s  research 
(Ret  o| 

B.  EXPERIMENTAL  DESIGN 

The  experiment  is  designed  to  run  the  (  S  model  for  ninety-nine  timesteps  over 
thirtv  replications,  to  generate  force  level  data  and  calculate  the  resulting  average 
attrition  rate  coefficients.  Ihese  average  coefficients  are  then  used  as  the  input 
parameters  to  the  M  W  model  which  then  provides  the  calculated  force  levels  per 
timestep  This  force  level  data  is  used  tor  analvsis  I  he  method  lor  analvsis  w  hich  has 
been  chosen  is  to  plot  the  resulting  force  level  data  over  time  to  provide  lor  an  easier 
comparison  than  sorting  through  the  enormous  amount  ol  data  generated  for  each  test 
case 

I  his  experiment  was  conducted  for  cases  involving  Ixl.  1x2.  2x2,  and  tx?  force 
level  scenarios,  [  samples  of  the  data  collected  and  the  resulting  force  level  trajectorv 
plots  are  provided  in  the  next  section 

i  Model  Verification 

Before  beginning  the  experiment  using  these  two  models,  each  ol  the  models 
was  exercised  and  tested  heavily  to  verifv  that  thev  ran  properlv  and  provided  the 
desired  output  as  they  were  designed  A  brief  description  ol  the  verification  process  tor 
each  model  will  be  given 

a.  Computer  Simulation  Model  K  /  S  model) 

I  he  first  step  in  veritving  the  <  S  model  involved  a  procedure  of  several 
runs  using  different  input  data  sets  to  test  the  program  to  ensure  that 


1)  The  input  data  and  program  were  compatible  and  the  input  variable*  are 
properly  read. 

2)  All  the  computational  code  was  performing  properly  (i.e..  the  average  attrition 
rate  coefficients  and  the  average  force  levels* 

3)  That  each  of  the  test  conditions  was  functioning  as  designed  n.e  .  the 
condition  for  ending  at  the  specified  force  levels  \SlOP  and  YSlOPl. 

4)  That  the  output  was  presented  as  desired  lor  the  step-h\-step  analvsis  and  the 
plotting  data. 

5 1  Test  the  resulting  average  attntion-rate  coefficients  with  the  input  probabilities 
of  lull  which  correspond  to  these  coefficients. 

During  this  portion  of  the  research  it  was  interesting  to  note  that  the  (  S 
model  was  designed  to  allow  each  shooter  one  shot  per  time  interval,  but  the  shooter 
could  select  a  target  to  lire  at  which  had  alrcads  been  lulled  during  that  time  interval 
I  his  fact  seemed  to  be  a  lairlv  realistic  condition  that  could  occur  on  the  battle  field 
where  two  combatants  lire  at  the  same  target  instantaneously  or  within  a  traction  ol  a 
second  ol  each  other  Therefore,  this  condition  was  left  in  the  model 
b.  Moose  j  Wozemraft  Model  <  Ml  W  model) 

A  computer  program  was  designed  to  run  the  M  W  model  to  allow  lor  an 
easy  comparison  of  the  two  models  Bs  utilizing  a  computer  program  it  allowed  lor 
the  same  timestep  by  timestep  analvsis  as  the  (  S  model 

After  designing  the  TORI  RAN  program,  shown  in  Appendix  D.  to  run  the 
Moose  Wozencrafl  model  iM  W  model!  it  was  tested  to  verily  that  the  following 
conditions  were  met: 

h  The  input  data  and  the  program  were  compatible  and  all  variables  were  read 
properly 

2)  All  computational  code  was  performing  as  designed. 

3>  The  resulting  output  files  provided  the  desired  data  lor  analysis  Since  the 
M  W  model  is  a  deterministic  model,  the  same  output  will  be  obtained  each 
time  for  the  same  set  of  input  data  parameters 
Each  of  the  conditions  listed  above  was  verified  for  simple  input  data  sets  for  each  of 
the  cases 

2  Data  Generation  Procedure 

The  procedure  followed  in  conducting  this  experiment  was  a  three  step  process 
lor  each  unique  case  tested.  I  he  first  step  involved  running  the  ('  S  model  for  ninetv- 
nme  timesteps  over  thirtv  replications  which  would  generate  the  following  data  files 

'I 


1)  The  calculated  average  attrition-rate  coefficients  which  are  later  used  as  input 
to  the  M  W  model.  An  example  file  is  shown  in  I  igure  3  I 
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A  VC  RAM  ATTRITION  COEFFICIENTS  FOR  Y  FORCE 

AAll.ll  AAIl.21  AAI1.SI  AAlf.ll  AAU.2I  AAlt.Sl  AAll.ll  AAU.2>  AAll.Sl 

090012  0.00000  0.00000  0  00001  0  00000  0  OOOOO  0  00000  0  00000  0  00000 

AVCRAM  ATTRITION  COCFFICICNTS  FOR  Y  FORCE 

Mil, II  Mil. 21  Mil. SI  MI2.1I  Ml  2  ,2  >  MI2.JI  Ml  1 . 1  l  MIS. 21  MIS, Si 

0  0*017  0.00000  0.00000  0 . 00142  0  OOOOO  0.00000  0.00000  0  OOOOO  0  OOOOO 

AVCRAOC  ATTRITION  COEFFICIENTS  FOR  X  FORCE 

ecu. it  ceil. 2i  ecu. si  ecu. it  ecu. 21  ccu.si  ccis.n  ecu. 21  ccu.si 

0  00027  0.00000  0.00000  0  OOOOO  0.00000  0  OOOOO  0  OOOOO  0  OOOOO  0  OOOOO 

AVCRAOC  ATTRITION  COEFFICIENTS  FOR  X  FORCE 

MI1, II  MU. 21  Mll. 31  MI2.ll  MU, 21  MU. Si  DOU.ll  MIS, 21  MU, Si 

0.00*40  0  OOOOO  0 . OOOOO  0  OOOOO  0  OOOOO  0  OOOOO  0  OOOOO  0  OOOOO  0  OOOOO 

STORPCO  RY  PROQRAM  AT  NBTOP  ■  M 

l  vgure  3  i  Lxample  Average  Ammon-rate  i  oefficients  l  tic  l  or  t  S  Model 

2l  A  lile  containing  the  average  force  level  lor  each  lorce.  the  average  force  level 
plus  one  standard  deviation,  and  the  average  force  level  minus  one  standard 

deviation  is  provided  for  the  aggregate  force  levels  \  t^x  l  and  Y  i^v  i  per 

timestep  as  well  as  for  each  of  the  individual  units  xf  and  v^  An  example  file 
is  shown  in  Figure  3  2 

3)  I  he  average  force  level  and  standard  deviation  of  each  force  level  per  time 
step.  An  example  of  these  results  is  shown  in  Figure  3.3 

The  second  step  involves  using  the  average  attrition-rate  coefficients  calculated 
by  the  C  S  model  as  input  parameters  to  the  M  W  model  I  he  output  generated  hv 
the  M  W  model  for  each  set  of  input  data  is  the  following 

1)  The  force  level  for  each  of  the  forces  x]  and  v^  as  well  as  the  total  force  levels 
X  <Vy  and  Y  (^y  (  lor  each  time  step  An  example  of  this  output  tile  is 
shown  in  Figure  3  4. 
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1  igure  3. 
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F  igure  3.3  Example  File  of  Force  Level  and  Standard  Deviation  1  or  (.'  S  Model 

The  third  step  involves  plotting  the  force  level  data  generated  by  each  of  the 
models  to  perform  the  comparison.  The  plots  which  were  utilized  for  the  comparison 
are  listed  below  with  examples  provided: 

1  >  A  plot  of  the  total  X  force  vs.  Total  Y  force  level  trajectories  for  both  models 
is  done  first.  This  plot  provides  for  a  macro  level  comparison  of  the  resulting 
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1  igurc  3  4  Example  \-*\i  1  orce  l  evel  Data  1  ile  1  or  \l  W  Model 

lorce  levels  over  time  lie,  the  lorce  (  ommanders  lev  cl  i  This  plot  is  two 
dimensional,  plotting  the  \  forte  level  vs  8  lorte  level  I  he  time  variable  is 
shown  bv  marking  the  lorte  level  everx  fifth  timestep  I  his  plot  will  be  used 
to  analwe  the  two  aggregate  forte  models  results  An  example  of  the  total 
forte  icvel  trajectory  plot  is  shown  in  1  igure  3  *> 

2i  Plots  ol  the  x(  lorte  level  vs  time  and  the  v ^  lorte  levels  vs  time  tor  both  the 
(  S  and  the  M  W  models  are  given  I  his  provides  a  comparison  ot  the  battle 
outcome  tor  the  two  models  at  the  micro  level  An  example  of  the  \)  tor  v^i 
plot  is  shown  in  T  igurc  3  6 

The  plots  introduced  in  this  chapter  are  used  extcnsivelv  tor  conducting  the 
comparison  ot  the  two  models  C  hapter  1  our  will  present  the  results  ot  the  experiment 
and  discuss  the  comparisons 


'4 


TOTAL  FORCE  LEVEL  TRAJECTORY 


Ml  MO 

400  4M  OOO  MO  000  MO  TOO 

TOTAL  X  FORCE  LEVEL 

CASE  10  (3X3) 

[tl  rT 


Tl 


XI  FORCE  LEVEL 

•  a  r\A 


8  O  i  10  i»  *'•  «»  s°  34  ,0  4i,u 


time 


CASE  10  (3X3) 


Figure  3  6 


,  0r  x  force  level  trajectories. 

.6  Example  plot  oi  Xj 


IV.  EXPERIMENTAL  RESULTS 


A.  INTRODUCTION 

The  primary  objective  of  this  experiment  was  to  test  whether  the  deterministic 
Moose;  Wozencraft  model  w'as  a  good  approximation  (or  prediction)  for  the  mean 
result  of  the  stochastic  computer  simulation.  The  secondary  objective  for  this 
experiment  w’as  to  test  whether  the  two  models'  force  level  trajectories  behaved  in  the 
expected  manner  described  by  Professors  Moose  and  Wozencraft.  The  results  of  ten 
cases  are  presented  graphically  showing  the  two  models'  results  with  a  bracket  of  plus 
and  minus  one  standard  deviation  as  a  measure  of  uncertainty.  Each  case  will  be 
discussed  in  the  context  of  the  primary  and  secondary  objectives  in  this  chapter  and 
their  results  shown  graphically  in  Appendix  E. 

B.  THE  FOUR  MODES  OF  BEHAVIOR 

Before  proceeding  with  the  presentation  of  the  results,  it  is  important  to  briefly 
explain  the  four  trajectory  behavior  modes  which  will  be  referenced  in  each  of  the 
cases.  Professors  Moose  and  Wozencraft  found  that  the  force  level  trajectories 
displayed  one  of  four  distinct  modes  of  behavior.  A  mode  of  behavior  was  the  result  of 
the  existence,  or  nonexistence,  of  equilibrium  points  in  the  state  space  (positive  first 
quadrant).  The  four  modes  related  to  the  existence  of  one  stable  equilibrium,  one 
unstable  equilibrium,  two  equilibria  (where  one  is  stable  and  one  is  unstable),  or  no 
equilibria.  These  four  modes  of  behavior  are  illustrated  in  Figure  4.1.  [Ref.  6] 

The  expected  behavior  for  any  force  level  trajectory  when  one  stable  equilibrium 
point  exists  in  the  state  space  is  shown  in  Figure  4.1  (a).  The  trajectory  will  be 
attracted  to  the  stable  equilibrium  and  remain  there,  since  dx/dt=  dy,'dt  =  0  at  the 
stable  equilibrium. 

The  expected  behavior  for  any  force  level  trajectory  when  one  unstable 
equilibrium  point  exists  in  the  state  space  is  shown  in  Figure  4.1  (b).  As  the 
trajectories  near  the  equilibrium  they  are  repelled  or  turned  away  from  the  equilibrium 
point.  The  presence  of  an  unstable  equilibrium  creates  a  division  of  the  state  space 
into  two  regions.  The  dividing  line  was  called  a  separatrix.  In  scenarios  which  exhibit 
this  mode  of  behavior  the  outcome  depends  on  which  side  of  the  separatrix  the  initial 
force  levels  begin. 
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The  expected  behavior  for  any  force  level  trajectory  when  two  equilibria  exist 
(one  stable  and  one  unstable)  is  shown  in  Figure  4.1  (c).  Once  again  the  unstable 
equilibrium  creates  a  separatrix  dividing  the  state  space  into  regions.  The  stable 
equilibrium  will  attract  trajectories  as  discussed  earlier.  This  particular  mode  of 
behavior  has  two  possible  outcomes  for  a  battle,  one  of  the  forces  can  win  or  they  will 
fight  to  a  draw.  The  outcome  depends  on  the  initial  force  levels  for  each  force. 

The  expected  behavior  for  any  force  level  trajectory  when  no  equilibrium  point 
exists  is  shown  in  Figure  4.1  (d).  The  trajectories  will  approach,  and  eventually  cross 
one  of  the  axes,  which  results  in  a  force  being  eliminated. 

C.  DISCUSSION 

This  numerical  experiment  was  designed  to  test  a  combat  scenario  involving 
different  opposing-force  composition.  (N„  x  Mv)  cases  involving  N'x  types  of  X  forces 
vs.  M,,  types  of  Y  forces  are  tested  in  this  experiment.  An  example  is  a  1x2  case  which 
involves  one  X  type  force  vs.  two  Y  type  forces.  The  input  data  sets  were  chosen  to 
test  the  two  models  for  the  four  distinct  modes  of  behavior  that  Professors  Moose  and 
Wozencraft  identified  [Ref.  6].  The  computed  results  are  presented  in  sections  which 
are  dedicated  to  a  specific  (X  x  M)  scenario. 

Appendix  E  contains  the  input  data  sets  used  for  each  model,  and  the  resulting 
force  level  trajectory  plots  for  each  of  the  cases.  In  each  case  the  input  data  for  the 
C,  S  model  was  developed  to  test  one  of  the  modes  of  behavior  and  the  M  W  model 
input  data  were  the  calculated  attrition  rate  coefficients  calculated  by  the  C/S  model. 
Several  trajectory  plots  were  done  for  each  of  the  cases.  There  is  a  plot  of  the  total 
aggregate  force  level  trajectory  which  compares  the  total  X  force  vs.  the  total  Y  force. 
This  aggregate  trajectory  shows  the  battle  outcome.  An  example  of  this  plot  is  shown 
for  case  six  in  Figure  4.2. 

The  result  of  the  computer  simulation  model,  which  is  the  mean  of  thirty 
replications,  is  always  shown  as  a  solid  line  with  squares  marking  the  force  level  every 
fifth  time  step.  The  Moose, 'Wozencraft  model  results  are  shown  as  a  dashed  line  with 
circles  marking  the  force  level  every  fifth  time  step.  Additionally,  two  solid  lines  are 
shown  wdth  no  markings.  These  are  the  curves  which  show  plus  and  minus  one 
standard  deviation  from  the  C/S  model  results.  They  serve  as  a  measure  of  uncertainty 
for  the  comparison.  The  aggregate  force  level  plot  provides  a  comparison  of  the  iwo 
models  total  force  level  results  to  determine  whether  they  do  have  the  same  general 
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behavior  and  whether  the  deterministic  model  is  a  good  approximation  for  the  result  of 
the  stochastic  model.  Each  of  the  aggregate  forces  (i.e.,  X  =  Yx-)  are  plotted  as  a 
function  of  time  to  provide  a  clearer  comparison  of  the  two  models  results  as  the  battle 
progresses.  An  example  of  these  plots  are  shown  in  Figures  4.3  and  4.4.  This  provides 
another  comparison  between  the  two  models'  results  for  each  aggregate  force.  In 
addition  to  the  aggregate  force  level  plots,  the  results  for  each  of  the  separate  units  are 
plotted  as  a  function  of  time.  Examples  of  these  plots  are  provided  in  Figures  4.5 
through  4.8.  This  provides  a  comparison  of  the  two  models  at  the  individual  unit  level. 

The  method  of  comparison  used  for  this  analysis  was  a  graphical  presentation  of 
the  force  level  trajectories  rather  than  an  extensive  numerical  analysis  of  the  resulting 
data.  Each  case  will  be  presented  in  a  similar  manner.  The  cases  will  be  briefly 
described  and  then  a  list  of  observations  given.  Specific  comments  are  made  about 
whether  the  M/W  trajectory  was  close  to  the  C/S  trajectory  and  whether  the  expected 
behavior  was  observed. 

D.  (1X1)  CASE  COMPARISONS 

The  experiment  was  conducted  for  cases  one  through  four  using  one  X  force  and 
one  Y  force,  which  are  labeled  as  lxl  cases.  The  results  of  each  of  these  cases  are 
shown  in  Figures  E.l  through  E.12  in  Appendix  E.  The  force  level  trajectories  for  the 
M/W  model  were  very  close  to  the  mean  force  level  trajectories  for  the  C.'S  model  in  all 
four  cases. 

1.  Case  One 

Case  one  was  designed  to  test  the  two  models  involving  a  lxl  scenario  for  the 
mode  of  behavior  where  no  equilibrium  exists.  The  following  observations  are  noted 
from  the  results  shown  in  Figures  E.l  through  E.3: 

a)  The  resulting  force  level  trajectories  for  the  M/W  model  and  the  C  S  model 
were  very  similar,  both  displaying  the  same  trajectory  shape. 

b)  The  M/W  model  trajectory  was  well  within  the  one  standard  deviation 
boundaries,  for  the  complete  battle.  In  fact,  the  two  curves  over  lapped 
during  portions  of  the  battle.  This  case  resulted  in  the  X  force  winning  the 
battle. 

2.  Case  Two 

Case  two  was  designed  to  test  the  mode  when  one  stable  equilibrium  exists  in 
the  vicinity  of  (200,200).  The  following  observations  are  noted  from  the  results  shown 
in  Figures  E.4  through  E.6: 


Figure  4.3  X  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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a)  The  M  W  trajectory  and  the  C  S  trajectory  were  very  similar  and  they  display 
the  same  trajectory  shape. 

b)  The  M  W  model  trajectory  was  outside  the  established  boundaries  for  the 
initial  portion  of  the  battle.  It  is  noted,  however,  that  the  area  bounded  by 
one  standard  deviation  is  extremely  small  in  this  region.  The  M  W  trajectory 
entered  the  bounded  area  as  the  two  curves  turned.  Both  models  reached  an 
equilibrium  point  in  the  vicinity  of  (200,200)  where  the  battle  would  end  in  a 
draw. 

c)  Examination  of  the  individual  unit  force  level  trajectories  in  Figures  E.5  and 
E.6  shows  that  the  M  W  trajectory  was  within  the  established  bounds  lor  each 
force  as  a  function  of  time. 

3.  Case  Three 

Case  three  was  designed  to  test  the  1x1  scenario  for  the  mode  of  behavior 
where  one  unstable  equilibrium  exists  in  the  vicinity  of  (200,200).  The  following 
observations  are  noted  from  the  results  shown  in  Figures  E.7  through  E.9: 

a)  The  force  level  trajectories  of  both  models  displayed  the  expected  behavior  for 
a  case  where  one  unstable  equilibrium  exists.  Each  of  the  models  trajectories 
approached  the  equilibrium  and  then  turn  away  from  it. 

b)  Once  again  the  M  W  trajectory  was  well  within  one  standard  deviation  of  the 
C  S  trajectory. 

4.  Case  Four 

Case  four  was  designed  to  test  a  lxl  scenario  where  the  mode  of  behavior 
involves  two  equilibria  in  the  vicinity  of  (200,200KL)  and  (366.100KS).  The  following 
observations  are  noted  from  the  results  of  this  case  shown  in  Figures  E.10  through 
E.  1 2: 

a)  The  results  show  the  trajectories  of  both  models  display  the  expected  behavior 
approaching  the  unstable  equilibrium  and  then  turning  away. 

b)  The  M/W  model  was  always  well  within  the  one  standard  deviation 
boundaries  for  this  case. 

E.  (1X2)  CASE  COMPARISONS 

One  case  was  developed  to  test  an  unbalanced  force  level  scenario  involving  one 
type  X  force  and  two  type  Y  force.  The  purpose  of  this  case  was  to  ensure  that  both 
models  perform  properly  for  an  unbalanced  number  of  forces  engaging  in  the  battle.  A 
1x2  scenario  was  chosen  for  this  case  and  the  results  are  presented  next. 


Case  five  was  the  one  unbalanced  scenario  developed  and  the  mode  of 
behavior  tested  involved  one  unstable  equilibrium  in  the  vicinity  of  (30U.200).  The 
following  observations  were  noted  from  the  results  shown  in  Figures  E.13  through 
E  18: 

a>  Both  models  displayed  the  expected  behavior  for  an  unstable  equilibrium.  The 
trajectories  for  both  models  approached  the  equilibrium  point  and  then  were 
turned  away. 

b)  The  \i  W  trajectory  was  well  within  the  one  standard  deviation  boundaries, 
and  had  the  same  shape  as  the  C  S  trajectory.  It  is  interesting  to  note  the 
relative  speed  at  which  the  two  trajectories  moved  away  from  the  unstable 
equilibrium  point.  The  C  S  model  moved  away  much  faster  than  the  M  W 
model. 

c)  The  plot  of  each  force  vs.  time,  shown  in  Figures  E.14  through  E.18.  show  the 
two  models  results  separating  over  time,  and  the  M  W  trajectory  leaving  the 
bounded  area  toward  the  end  of  the  battle. 

d)  This  case  involving  a  1x2  scenario  with  an  unstable  equilibrium  is  similar  to 
case  three  which  was  a  lxl  scenario  involving  one  unstable  equilibrium. 

F.  (2X2)  CASE  COMPARISONS 

The  next  step  in  the  experiment  was  to  develop  cases  which  involved  two  types  of 
forces  on  each  side  and  compare  the  models  results  for  the  four  modes  of  behavior 
discussed.  Cases  six  through  nine  are  2x2  case  scenarios.  Their  results  are  shown  in 
Figures  E.19  through  E.46.  The  results  for  each  of  these  cases  will  be  discussed 
individually. 

1.  Case  Six 

Case  six  is  a  2x2  force  level  scenario  designed  to  test  the  models  when  there 
are  no  equilibria  in  the  state  space.  The  results  for  case  six  are  shown  in  Figures  E.19 
through  E.25.  The  following  observations  are  noted  from  the  results: 

a)  The  resulting  trajectories  shown  in  the  Total  Force  Level  Trajectory  plot  in 
Figure  E.19  show  both  models  behaved  similarly.  In  this  particular  case  the  V 
force  was  the  winning  force.  The  expected  behavior  for  no  equilibria  was 
observed  in  this  case. 


b)  The  M;W  model  results  were  within  the  defined  one  standard  deviation 
boundaries  during  most  of  the  battle,  leaving  the  boundaries  only  toward  the 
end  of  the  engagement. 

c)  Further  comparison  of  the  two  models  can  be  done  using  the  results  plotted 
for  the  individual  forces  in  Figures  E.20  through  E.25.  From  these  results  it  is 
evident  that  the  X  force  results  are  the  ones  where  the  models  differed  the 
most.  The  trajectories  all  had  the  same  general  trends. 

d)  This  2x2  case  involving  no  equilibria  is  similar  to  case  one  which  was  a  1x1 
case  scenario  involving  no  equilibria. 

2.  Case  Seven 

Case  seven  was  designed  to  test  the  models  results  with  two  type  X  forces  and 
two  type  Y  forces  when  two  equilibria  exist.  The  results  of  the  two  models  are  shown 
in  Figures  E.26  through  E.32.  The  following  observations  are  noted  for  case  seven: 

a)  The  resulting  trajectories  for  both  models  behaved  in  the  expected  manner. 
They  approached  the  unstable  equilibrium  and  then  turned  away  and  headed 
off.  In  this  particular  case  the  X  force  is  winning,  but  the  results  would  be 
different  if  the  force  levels  were  such  that  the  initial  point  was  above  the 
separatrix. 

b)  The  M/W  trajectory  had  the  same  shape  as  the  C/S  trajectory.  The  M/W 
trajectory  showm  in  Figure  E.26  was  outside  the  established  boundaries  during 
the  initial  portion  of  the  battle,  but  wras  within  the  boundaries  after  they 
turned  away  from  the  equilibrium.  However,  examination  of  the  individual 
units  force  level  trajectories  shown  in  Figures  E.27  through  E.32  shows  the 
M/W  model  within  the  boundaries  for  the  X  force  and  just  outside  the 
boundaries  for  the  Y  force. 

c)  This  2x2  case  involving  two  equilibria  is  similar  to  case  four  which  was  a  lxl 
scenario  involving  two  equilibria. 

3.  Case  Eight 

Case  eight  was  designed  to  test  the  two  models  results  for  a  force  scenario 
involving  two  type  X  forces  and  two  type  Y  forces  with  one  stable  equilibrium  present 
in  the  state  space.  The  results  for  case  eight  are  shown  in  Figure  E.33  through  E.39. 
The  following  observations  are  noted  for  case  eight: 

a)  The  resulting  trajectories  for  the  M,W  and  C/S  models  both  exhibited  the 
expected  behavior  with  one  stable  equilibrium  present.  They  both  were 
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attracted  to  the  stable  equilibrium  point  and  remained  in  its  vicinity.  The  C  S 
model  shows  more  variability  due  to  its  stochastic  process. 

b)  The  M/W  trajectory  had  the  same  shape  and  remained  within  the  established 
boundaries  for  the  entire  battle. 

c)  This  2x2  case  involving  one  stable  equilibrium  is  similar  to  case  two  which  was 
a  lxl  scenario  involving  one  stable  equilibrium. 

4.  Case  Nine 

Case  nine  was  designed  to  test  the  two  models  results  for  a  2x2  force  level 
scenario  when  one  unstable  equilibrium  exists.  The  results  for  case  nine  are  shown  in 
Figures  E.40  through  E.46.  The  following  observations  are  noted  for  case  nine: 

a)  The  resulting  trajectories  for  each  of  the  models  behaved  in  the  expected 
manner  for  an  unstable  equilibrium. 

b)  The  M/W  trajectory  had  the  same  general  shape,  approaching  and  turning 
away  from  the  equilibrium,  but  it  was  not  within  the  designated  boundaries 
during  most  of  the  battle.  Reviewing  the  results  of  the  individual  units  shows 
that  the  tw  o  models  did  not  behave  similarly  as  a  function  of  time. 

c)  This  2x2  case  involving  one  unstable  equilibrium  is  similar  to  case  three  which 
is  a  lxl  scenario  involving  one  unstable  equilibrium. 


G.  (3X3)  CASE  COMPARISONS 

The  last  force  level  scenario  developed  for  this  experiment  involves  three  type  X 
forces  and  three  type  Y  forces.  One  case,  labeled  as  a  3x3  case,  was  developed  for  this 
experiment.  There  is  no  discussion  of  the  two  models  results  for  each  of  the  four 
modes  of  behavior  discussed  earlier  for  the  3x3  scenario.  Rather,  one  case  was  run  to 
test  the  models  results  for  the  larger  3x3  scenario. 

1.  Case  Ten 

Case  ten  is  the  one  3x3  case  developed  for  comparing  the  M/W  model  results 
with  the  C/S  model  results.  The  following  observations  are  made  from  the  results 
showm  in  Figures  E.47  through  E.55: 

a)  The  M/W  trajectory  was  well  within  the  established  boundaries  throughout 
the  battle  and  the  two  trajectories  shapes  were  almost  identical. 

b)  The  M/W  and  C/S  trajectories  for  each  of  the  individual  units  are  very  similar. 
In  several  of  the  cases  the  two  trajectories  are  overlapping  for  most  of  the 
battle. 
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H.  SUMMARY  OF  RESULTS 

The  results  of  the  numerical  experiment,  with  ten  cases  involving  lxl,  1x2,  2x2, 
and  3x3  force  level  scenarios,  have  been  presented  and  briefly  discussed  in  this  chapter. 
Some  general  observations  from  these  results  will  serve  as  a  summary  of  the 
experiment. 

First,  the  results  for  the  ten  cases  support  the  idea  that  a  deterministic  attrition 
model  is  often  a  good  approximation  for  the  mean  result  of  a  stochastic  attrition 
model.  From  the  results  shown  in  the  previous  section,  only  one  of  the  ten  cases 
displayed  force  level  trajectories  where  the  M/W  model  was  significantly  different  from 
the  C/S  model.  Five  of  the  ten  cases  showed  the  M/W  trajectories  were  within  the 
established  one  standard  deviation  boundary  for  the  entire  battle  (cases  one,  three, 
four,  eight,  and  ten).  Of  the  remaining  cases  where  the  deterministic  M  W  model 
results  were  outside  the  boundary  for  a  portion  of  the  battle  shown  in  the  Total  Force 
Level  Trajectory  plot,  it  was  noted  that  the  individual  unit  force  level  trajectories  were 
within  the  established  one  standard  deviation  boundary  in  all  but  one  case  (case  nine). 

Secondly,  the  results  from  the  experiment  show  that  both  models  force  level 
trajectories  displayed  the  expected  mode  of  behavior  for  the  number  and  type  of 
equilibrium  points  which  were  present  in  all  cases.  This  shows  strong  support  for  the 
research  done  by  Professors  Moose  and  Wozencraft. 

This  experiments  results  are  certainly  not  proof  that  the  deterministic  model  is 
always  a  good  approximation  for  the  mean  result  of  a  stochastic  model.  However,  the 
results  do  indicate  they  often  may  be  a  good  approximation. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  results  of  this  research  showed  that  the  M/W  trajectory  was  a  good 
approximation  for  the  C/S  trajectory  in  nine  of  ten  cases.  This  is  an  indication  that 
deterministic  attrition  models  may  be  good  approximations  for  the  mean  of  stochastic 
attrition  models.  It  has  also  shown  that  the  modes  of  behavior  for  both  models 
trajectories  are  very  similar.  The  combination  of  these  two  results  provides  strong 
support  that  the  M/W  model  is  a  good  approximation  for  the  C  S  model.  These 
results  provide  important  support  for  the  argument  that  a  deterministic  attrition  model 
can  be  used  to  model  the  combat  attrition  process  which  is  generally  agreed  upon  as  a 
stochastic  process. 

This  simulation  provided  a  test  of  the  Moose/Wozencraft  model.  The  M/W 
model  was  tested  utilizing  four  different  force  composition  scenarios,  for  each  of  the 
four  types  of  trajectory  behavior  modes.  The  results  showed  the  M/W  models  results 
were  good  approximations  for  the  individual  unit  level  and  the  aggregate  force  levels  of 
the  C/S  model.  This  Generalized  Lanchester  Equation  type  model  is  simple  to  use  and 
understand.  It  allows  the  modeler,  or  military  planner,  to  model  an  aggregate  force  of 
as  many  different  type  forces  as  required  to  model  the  real-world  battlefield.  The  M/W 
model  has  real  potential  for  future  use  by  military  analysts. 

B.  COMMAND  AND  CONTROL  APPLICATIONS 

A  key  decision  in  the  Command  and  Control  decision  process  is  what  force 
structure  is  required  to  engage  a  proposed  enemy.  The  problem  is  to  determine  the 
optimum  force  size  and  composition  required  to  engage  an  enemy  and  win.  Optimizing 
one's  force  is  a  very  complex  problem,  involving  the  right  mixture  of  force  types  and 
the  right  number  of  each  of  the  forces.  An  example  would  be  when  a  commander  may 
need  a  tank  heavy  unit  for  a  specific  battle  in  order  to  ensure  that  his  force  will  win. 

A  model  such  as  the  M/W  model  can  serve  as  a  very  useful  decision  aid  for  the 
military  planners  and  combat  commanders  to  assist  with  optimization  of  forces  on  the 
battlefield.  Given  an  accepted  set  of  attrition-rate  coefficients,  the  model  can  provide 
the  capability  of  testing  for  the  optimum  mix  and  force  size  required  to  fight  and  win  a 
battle.  The  use  of  a  model  of  this  type  by  maneuver  units  would  require  the 
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commanders  and  their  stafTs  to  ask  the  intelligence  community  for  specific  information 
on  the  enemies  current  force  structure  and  resupply  capability.  This  information  could 
be  used  with  the  model  to  assist  planners  and  commanders  to  decide  whether  to 
engage,  and  the  optimum  force  compostion  if  the  decision  is  made  to  do  so.  If  not, 
they  may  chose  to  with-draw  and  fight  on  better  terms.  Models  of  this  type  would  be 
beneficial  to  military  planners  responsible  for  developing  OPLANS  and  CONPLANS 
also.  It  would  allow  them  to  hypothesize  the  expected  enemy  force  structures  and  test 
proposed  force  structures  to  engage  the  enemy. 

C.  RECOMMENDATIONS 

After  completing  this  project  it  is  apparent  that  a  great  deal  of  future  research 
can  be  done  with  these  two  models  to  provide  a  more  complete  analysis.  The  results 
from  this  initial  survey  have  identified  areas  for  further  testing  which  would  provide  a 
more  complete  understanding  of  the  dynamics  of  the  two  models.  The  areas  listed 
below  would  provide  good  topics  for  future  work  with  these  models: 

1)  Utilize  the  models  to  locate  the  equilibrium  points  and  the  separatrix  for  cases 
where  an  unstable  equilibrium  point  exits  in  the  state  space.  Once  they  have 
been  located,  observe  the  tw'o  models  results  in  the  vicinity  of  the  equilibrium 
point  for  several  different  inital  force  levels.  The  separatrix  is  a  very- 
interesting  factor  which  should  be  observed  closely. 

2)  A  detailed  analysis  of  case  nine  of  this  experiment  is  recommended.  Case  nine 
was  the  one  case  which  was  run  that  did  not  behave  similarly  to  the  rest  of  the 
cases.  At  this  time  it  is  not  known  why,  and  therefore  deserves  further 
research. 

3)  One  of  the  observations  made  from  this  experiment  involved  a  difference  in 
the  attrition-rate  of  a  force  over  time  for  the  two  models.  This  could  be  a 
result  of  the  normalization  factor  used  for  the  aimed  fire  attrition  in  the  M  W 
model.  This  factor  was  used  to  ensure  that  an  individual  force  was  not 
allowed  to  fall  below  a  zero  level.  This  factor  may  cause  the  two  models  to 
have  slightly  different  attrition-rates  over  time.  This  can  be  tested  by- 
removing  the  factor  and  changing  the  models  computer  program  to  stop  if  one 
force  falls  to  a  zero  level. 

4)  Test  the  model(s)  against  approved  combat  simulations. 
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Further  research  in  this  area  may  lead  to  an  attrition  model  which  is  accepted  as 
a  useful  model  for  modem  combat.  This  could  provide  a  useful  Command  and  Control 
decision  aid  for  use  by  commanders  and  their  staffs  to  optimize  force  structures  on  the 
battlefield. 


APPENDIX  A 

EXEC  FILE  FOR  C/S  PROGRAM 


&TRACE  ON 

GLOBAL  TXTLIB  VFORTLIB  NONIMSL 
GLOBAL  LOADLIB  VFLODLIB 
FILEDEF  4  TERM 
FILEDEF  5  DISK  INTAW  DATA  A 
FILEDEF  6  DISK  TAW  FT06F001  A 
FILEDEF  7  DISK  TAW  FT07F001  A 
FILEDEF  8  DISK  TXTY  DATA  A 
FILEDEF  9  DISK  X1Y1  DATA  A 
FILEDEF  10  DISK  X2Y2  DATA  A 
FILEDEF  11  DISK  X3Y3  DATA  A 
FILEDEF  12  DISK  STEP  DATA  A 
FILEDEF  13  DISK  STATTF  DATA  A 
FILEDEF  14  DISK  STATF1  DATA  A 
FILEDEF  15  DISK  STATF2  DATA  A 
FILEDEF  16  DISK  STATF3  DATA  A 
LOAD  TAW 
START  MAIN 

CP  MSG  0731P  JOB  TAW  COMPLETED 


IMSLSP  CMSLIB 
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APPENDIX  B 

C/S  COMPUTER  PROGRAM 


PROGRAM  LANCHESTER 


6/18/86 


WRITTEN  TO  TEST  THE  MOOSE/WOZENCRAFT  THEORY 


NTYPX , Y 
NXS , YS ( I ) 
NX , Y ( I ) 


NUMBER  OF  TYPES  OF  COMBATANTS  <INPUT> 
INITIAL  NUMBER  OF  EACH  TYPE  OF  COMBATANT 
NUMBER  OF  EACH  TYPE  OF  COMBATANT  ALIVE 


C 

AAA 

I,J 

ATTRITION: 

X 

'1 

c 

BBB 

l.J 

ATTRITION: 

Y 

J 

c 

CCC 

1  > J 

ATTRITION: 

Y 

I 

c 

DDD 

I,J 

ATTRITION: 

X 

J) 

LOSSES  TO  AREA  FIRE 
SHOOTS  X 

LOSSES  TO  AREA  FIRE 
SHOOTS  Y 


NRRR(I) 

NSSS(I) 


NUMBER  COMBATANT  X( 
NUMBER  COMBATANT  Y( 


RESUPPLIED/TIMESTEP 

RESUPPLIED/TIMESTEP 


X, YPOS (K) 

X,YSTAT(K) 

NQX,Y(K) 


STORED  POSITION  OF  EACH  COMBATANT 

STORED  STATUS  (LIVE  OR  DEAD)  OF  EACH  COMBATANT 

STORED  TYPE  OF  EACH  COMBATANT 


PKHX ,  Y  ( I ,  J ) 
X,YAF(l , J) 


PK:  X  SHOOTS  Y  TYPE  (I,J)  INTERACTION  <INPUT> 
AREA  FIRE  LOSS  COEFFICIENTS  FOR  X  <INPUT> 


RN(K) 

RX , Y ( K ) 

RXS ,  YS  ( K) 

RXT, YT(K) 


LISTS  OF  RANDOM  VARIANTS,  CONTROLS  WHO  SHOOTS 
LISTS  OF  RANDOM  VARIANTS,  CONTROLS  HITS 
LISTS  OF  RANDOM  VARIANTS,  CONTROLS  SELF-LOSSES 
LISTS  OF  RANDOM  VARIANTS,  CONTROLS  TARGET  CHOICE 


PKX,Y(I,J) 
KAX,Y(I,J 
KHX,Y(I , J ) 
PKSX,Y( I ) 


PK:  X(I)  KILLS  Y(J)  THIS  TIMESTEP 
NR  OF  Y(  J )  HIT  BY  X(I)  THIS  TIMESTEP,  AREA  FIRE 
NR  OF  Y( J)  HIT  BY  X(l)  THIS  TIMESTEP,  AIMED  FIRE 
PROBABILITY  OF  SELF-LOSS 


NDHX ,  Y  ( I ' 
NDAX ,  Y  ( I 

ndsx,y(i; 


X,Y(I)  LOSSES  TO  AIMED  FIRE  THIS  TIMESTEP 
X, Y( I )  LOSSES  TO  AREA  FIRE  THIS  TIMESTEP 
X,Y(I)  SELF  LOSSES  THIS  TIMESTEP 


NSTOP 

X,YLOC 


MAXIMUM  NUMBER  OF  TIMESTEPS  <INPUT> 
POSSIBLE  POSITIONS  OF  COMBATANTS  <INPUT> 


IMPLICIT  REAL*8  (D) 


CHARACTER*40  I HE AD 
CHARACTER*!  IANS,IYES 


DIMENSION  XPOS(32^2) ,XSTAT(3222) , NQX(3222) ,NX(3) 
DIMENSION  YPOS (3222) , YSTAT(3222) ,NQY(3222) ,NY(3) 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


XLOC(; 
YLOCh 
XAF(3, 
PKX(3, 
PKHX( • 

pkax(: 


,PKSX(3) ,NRRR(3) 
,PKSY(3) ,NSSS(3) 


,YAF  3,1 


i) , PKHY (3,3), KHX (3,3), KHY (3,3) 
3 }  ,  P  KA  Y  ( 3 , 3 ) ,  KA  X  ( 3 , 3 ) ,  KA  Y  ( 3 , 3 ) 
NDY ( 3 ) , NXS (3) , NYS ( 3 ) 


,NDHY(3) 

'  , NDAY ( 3 ) 

'  ,NDSY(3) 

I) , AAAAV (3,3), SAAAAV (3,3) 


non  on 


c 

c 


DIMENSION  BBB(3,3) ,BBBAV(3 , 3) , SBBBAV(3 , 3) 

DIMENSION  CCC(3 , 3 ) , CCCAV (3,3), SCCCAV (3,3) 

DIMENSION  DDD ( 3 , 3 ) , DDDAV (3,3), SDDDAV (3,3) 

DIMENSION  RX(4222) .RY(4222) ,RN(9888) 

DIMENSION  RXS(4222) ,RYS(4222) 

DIMENSION  RXT(4222) ,RYT(4222) 

DIMENSION  SNTX(212) ,SNX1(212) ,SNX2(212) ,SNX3(2I2) 
DIMENSION  SNTY(212),SNYl(212) ,SNY2(212) ,SNY3(212) 
DIMENSION  SSNTX(212) ,SSNX1(212) ,SSNX2(212) ,SSNX3(212) 
DIMENSION  SSNTY(212) ,SSNYl(212) ,SSNY2(212) ,SSNY3(212) 
DIMENSION  ANTX(212) , ANX1 (212) , ANX2 (21Z) , ANX3 (212) 
DIMENSION  ANTY(212) .  ANY1 (212 ) , ANY2 ( 212 ) , ANY3 (212 ) 
DIMENSION  SDNTX ( 212 ) , SDNX1 ( 2 1 2 ) , SDNX2 ( 2 12), SDNX3  212) 
DIMENSION  SDNTY(212 ) , SDNY1 (212 ) , SDNY2 (212 ) , SDNY3 (212) 
DIMENSION  PSNTX(212) ,PSNXl(212) ,PSNX2(212) ,PSNX3(212) 
DIMENSION  PSNTY(212) ,PSNYl(212) ,PSNY2(212  ,PSNY3(212 
DIMENSION  TSNTX(212) ,TSNXl(212) ,TSNX2(212  ,TSNX3(212) 
DIMENSION  TSNTY(212) ,TSNYl(212) ,TSNY2(212) ,TSNY3(212) 
DIMENSION  CSTEP(lOO) 

DATA  IYES/'Y'/ 


CALL  EXCMS ( ' CLRSCRN ' ) 

READ  (5,9400)  I HE AD 
WRITE (4 , 9410 ) 

9410  FORMATS  DO  YOU  WANT  TO  CHANGE  THIS  HEADING?  (Y/N)  (A40) 1 ) 
WRITE (4, 9400)  I HEAD 
READ  (4,9405)  IANS 
9405  FORMAT (Al) 

IF(IANS.EQ.IYES)  READ  (4,9400)  I HEAD 

9400  FORMAT (A40) 

ASK  USER  FOR  THE  NUMBER  OF  REPLICATIONS  DESIRED 
WRITE (4, 9490) 

9490  FORMATS  HOW  MANY  REPLICAS  DO  DESIRE?  (MAX  100)') 

READ(4 , 9495 )  L 
9495  FORMAT (14) 

WRITE (4, 9480) 

9480  FORMAT^  DO  YOU  WANT  FILES  6  AND  7  PRINTED?  ') 

READ  (4,9401)  IANS 

9401  FORMAT (Al) 

WRITE (6, 9400)  I HEAD 
WRITE (7. 9400)  IHEAD 
WRITE (8, 9400)  IHEAD 
WRITE (9, 9400)  IHEAD 
WRITE (10, 9400)  IHEAD 
WRITE (ll, 9400)  IHEAD 
WRITE(l2,9400)  IHEAD 
WRITE (13, 9400)  IHEAD 
WRITE ( 14, 9400)  IHEAD 
WRITE ( 15, 9400)  IHEAD 
WRITl(l6,9400)  IHEAD 
READ  (5,9500)  NSTOP 

READ  (5,9500)  NTYPX,NX,XLOC,PKSX,NRRR 
IF(IANS . EO . IYES )THEN 
WRITE(6 , 9*02) 

WRITE ( 6 , 9  50 1 )  NTYPX , NX , XLOC , PKSX , NRRR 
END  IF 

READ  (5,9505)  ( (PKHX( I ,J) , J=1 , 3) ,1=1,3) 

IF ( IANS. EO. IYES )THEN 
WRITE ( 6 , 9521 ) 

WRITE (6, 9520)  ( ( PKHX( I,J),J=1,3),I=1,3) 

ENDIF 

READ  (5,9505)  ( (XAF ( I , J ) , J=1 , 3 ) , 1=1 , 3 ) 

IF( IANS . EQ . IYES )THEN 


WRITE(6,9520)  ( (XAF(I , J) , J=1 , 3 ) , 1=1 , 3) 

END  IF 

READ  (5,9500)  NTYPY,NY,YLOC,PKSY,NSSS 

IF ( IANS. EQ. IYES) THEN 

WRITE(6,9503) 

WRITE(6,9501)  NTYPY,NY,YLOC,PKSY,NSSS 
ENDIF 

READ  (5,9505)  ( (PKHY(I, J) , J*1 ,3) , 1=1 ,3) 

IF ( IANS . EQ . I YES ) THEN 
WRITE (6 , 9522) 

WRITE(6 ,9520)  ( (PKHY( I , J) , J*1 , 3) , 1=1 , 3) 

ENDIF 

READ  (5,9505)  ( (YAF(I , J) , J=I ,3 ) , 1=1 , 3) 

IF ( IANS . EQ . IYES )THEN 
WRITE(6 , 9524) 

WRITE(6,9520)  ( (YAF(I , J) , J=l,3) ,1=1 ,3) 

ENDIF 

FORMAT (41 5 , 3F5.1,3F7.3,3I5) 

FORMAT(4I5 . 3F5 . 1 , 3F7 . 3 , 315 ) 

FORMAT (  '  #  X  TYPES  NUMBER  EACH  LOCATIONS 
X  1  RE -SUPPLY  1 ) 


9503  FORMAT (  '  #  Y  TYPES  NUMBER  EACH 
X  1  RE -SUPPLY  1 ) 

9505  FORMAT(3F8. 5) 

9520  FORMAT (9F8. 5,/) 

9521  FORMAT (/ , 1  KILL  PROBABILITIES  FOF 

X  1  PK (1 , 1 )  PK(1 ,2)  PK(1, 

X  '  PK(3 , 1)  PK(3 , 2)  PK(3 , 

9522  FORMAT (/ , '  KILL  PROBABILITIES  FOF 

X  1  PK( 1 , 1 )  PK(1 ,2)  PK(1 , 

X  1  PK(3 , 1 )  PK(3,2)  PK(3, 

9524  FORMAT (/ , 1  AREA  FIRE  COEFFICIENTS 
X  1  AA(l.l)  AA(1 , 2 )  AA( 1 , 

X  1  AA(3 , 1 )  AA ( 3 , 2 )  AA(3, 

9523  FORMAT (/ , 1  AREA  FIRE  COEFFICIENTS 

X  '  CC(1,1)  CC(1 , 2)  CC(1, 

X  '  CC(3 , 1 )  CC( 3,2)  CC(3, 


LOCATIONS 


PK(2,1)  PK(2,2) 
') 

ORCE  Y( I )  SHOOTS 
PK(2,1)  PK(2,2) 
') 


AA(2 , 1 )  AA(2,2)  AA 

DRCE  X(I)  SHOOTS  Y( 
CC(2,1)  CC(2 ,2)  CC 


SELF-LOSS 

SELF-LOSS 

f  WIW: 

I  kiU,  3)'\  ' 


INITIALIZE  THE  RANDOM  NUMBER  GENERATORS 

READ ( 5 , 9510)  DX , DY , DN 
DSEEDX=DX 
DSEEDY=DY 
DSEEDN=DN 

I F ( IANS . E0 • IYES ) THEN 
WRITE(6,9511) 

WRITE ( 6 , 9  5 10 )  DX , DY , DN , DSEEDX , DSEEDY , DSEEDN 
ENDIF 

FORMAT (3D15 .8 , 3111) 

FORMAT ( 1  DOUBLE  PRECISION  RANDOM  VARIANT  SEEDS 
X  •  INTEGER  SEED  VALUES 1 ) 


INITILAIZE  THE  STATISTICS  VARIABLES 

DO  10  1*1,  212 
SNTX( I )=0 
SNX1(I)=0 
SNX2(I)=0 
SNX3(I)=0 
SNTY(I)=0 
SNY1 (I )=0 
SNY2(I)=0 
SNY3 ( I )=0 
SSNTX( I )=0 
SSNX1  I  =0 
SSNX2 ( I ) =0 
SSNX3  I)=0 
SSNTY  I)=0 


non*-* 
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SSNY1 ( I )=0 
SSNY2 ( I  )=0 
SSNY3 ( I  )=0 
ANTX(I)=0 
ANX1 ( I  )=0 
ANX2(I)=0 
ANX3 (I )=0 
ANTY ( I  )=0 
ANY1 ( I )  =0 
ANY2(I)=0 
ANY3(I)=0 
SDNTX(I)=0 
SDNX1  I  =0 
SDNX2(I  =0 
SDNX3(I  =0 
SDNTY(l)=0 
SDNY1 ( I )=0 
SDNY2 ( I )  =0 
SDNY3 ( I  )=0 
PSNTXi I )=0 
PSNX1 ( I  j=0 
PSNX2(l)=0 
PSNX3(I)=0 
PSNTY ( I  )=0 
PSNY1 ( I  j=0 
PSNY2 ( I  )=0 
PSNY3(I)=0 
TSNTXiI )=0 
TSNX1 ( I )=0 
TSNX2(I)=0 
TSNX3(l)=0 
TSNTY ( I  )=0 
TSNY1 ( I  )=0 
TSNY2 ( I  }=0 
TSNY3(I)=0 
CONTINUE 
DO  15  1=1,3 
NYS?I)=NY(I) 
NXS(I)=NX(I) 

15  CONTINUE 
MINST=99 
AVSTEP=0 . 0 
SDSTEP=0 . 0 
SSTEP=0 . 0 
SSSTEP=0 . 0 
DO  17  1=1,100 
CSTEP(I)=0.0 
CONTINUE 

BEGIN  REPLICATIONS 

DO  4000  N=1,L 
NSTEP*0 
NTX=0 
NTY=0 

DO  20  1=1,3 


PX)  THEN 
NTX=NTX+NX( I ) 

END  IF 

IF(I.LE.NTYPY)  THEN 
NTY=NTY+NY ( I ) 

END  IF 

20  CONTINUE 
NTXS=NTX 
NTYS=NTY 
NXSTOP=50 
NYST0P=50 


v.y-w, 
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INITIALIZE  THE  TWO  FORCES 

IMAX=0 

IMIN=1 

IP=0 

DO  200  J=1 ,NTYPX 
IMAX=IMAX+NX( J) 

DO  100  I=IMIN, IMAX 
IP=IP+1 

XPOS ( IP)=XL0C( J) 

XSTAT ( IP)=1 . 0 
NQX(IP)=J 
CONTINUE 
IMIN=IMAX+1 
CONTINUE 

IKAX=0 

IMIN=1 

IP=0 

DO  400  J=1,NTYPY 
IMAX=IMAX+NY(J) 

DO  300  I=IMIN, IMAX 
IP=IP+1 

YPOS( IP)=YLOC( J) 

YSTAT(IP)=1.0 
NQY(IP)=J 
CONTINUE 
IMIN=IMAX+1 
CONTINUE 
DO  500  1=1,3 
NDX( I )=0 
NDY(I)=0 
CONTINUE 
DO  700  1=1,3 
DO  700  J=1 , 3 
AAAAV(I , J)=0 .0 
BBBAV< I ,J )=0.0 
CCCAV ( I , J)=0 . 0 
DDDAV ( I , J)=0 . 0 
CONTINUE 

SETUP  A  FILE  FOR  FUTURE  PLOTTING 

IF ( IANS. EO.IYES) THEN 
WRITE(7 ,9/00) 

FORMAT (•  #  NTXS  NTX  NX1  NX 2  NX3  DX1  DX2  DX3 ' 

X  '  NTYS  NTY  NY1  NY2  NY3  DY1  DY2  DY3 1 ) 

WRITE ( 7 , 97 1 0 )  NS TEP , NTXS , NTX , NX , NDX , NTYS , NTY , NY , NDY 
FORMAT ('+'14,2(515,314),/) 

WRITE (6  9711) 

FORMAT^#  NTXS  NTX  NX1  NX2  NX 3  RX1  RX2  RX3 1 , 

X  1  NTYS  NTY  NY1  NY2  NY3  RY1  RY2  RY3 1 ) 

WRITE (6 ,9710)  NSTEP , NTXS , NTX , NX , NRRR , NTYS , NTY , NY , NSSS 
END  IF 

NTT  TOTAL  NUMBER  OF  LIVE  COMBATANTS 

NSX.Y  NUMBER  OF  COMBATANT  WHO  IS  FIRING 

IX, Y  TYPE  OF  COMBATANT  FIRING 

NTARX, Y  COUNT  OF  TARGETS  ATTACKED 

NCX,Y  NUMBER  OF  TARGET  UNDER  FIRE 

JX , Y  TYPE  OF  TARGET  UNDER  FIRE 

THIS  IS  THE  BEGINNING  OF  THE  COMBAT  CYCLE 

CONTINUE 
NSTEP=NSTEP+1 
NTT=NTX+NTY 
NGT=1 . 5*REAL(NTT ) 


v.v.n  /- 


nnnnnnnnnnnnnnnnnnn 


NGX=1.3*REAL(NTX) 

NGY= 1 . 3  *REAL ( NTY ) 

C 

C  MAKE  THE  LISTS  OF  RANDOM  NUMBERS 

C 

CALL  GGUBS (DSEEDN , NGT , RN ) 

CALL  GGUBS (DSEEDX, NGX,RX) 

CALL  GGUBS (DSEEDX ,NGX ,RXS) 

CALL  GGUBS ( DSEEDX, NGX,RXT) 

CALL  GGUBS (DSEEDY,NGY,RY) 

CALL  GGUBS (DSEEDY,NGY,RYS) 

CALL  GGUBS ( DSEEDY , NGY , RYT ) 

C  ZERO  STORAGE  FOR  THIS  TIMESTEP 

C 

NSX=0 

NTARX=0 

NSY=0 


NTARY=0 
KSTOPX=0 
KSTOPY=0 
DO  1040  1=1,3 
NDX(I)=0 
NDY ( I  )=0 
NDHX( I )=0 
NDHY< 

NDAX( 

NDAY( 

NDSX( 

NDSY  ( 

DO  1( 

KAX( 

KAY 
KHX 

KHY(_ , 

1020  CONTINUE 
1040  CONTINUE 

DO  1060  1=1,3 
DO  1060  J=1 , 3 
PKAX(I,JJ=XAF(I,J] 


PKAY ( I , J 
1060  CONTINUE 
1080  CONTINUE 


:YAF(I,j; 


*NTY 

*NTX 


DO  1090  J=1 , 3 
DO  1090  1=1,3 
PKX(I , J)=PKAX 
PKY(I, J)=PKAY 
1090  CONTINUE 


+PKHX( I , J) 
+PKHY ( I , J ) 


DO  AIMED  AND  AREA  FIRE  INTERACTIONS  FOR  THIS  TIMESTEP 

CHOOSE  THE  COMBATANT  WHO  SHOOTS,  USING  RN(I) 

MAKE  COMBATANT  NUMBER  (COUNT),  NSX,Y 
COMBATANT  MUST  BE  ALIVE 
DETERMINE  COMBATANT  TYPE,  LX,Y 
COMBATANT  MUST  NOT  BE  A  SELF-LOSS 


INCREASE  THE  TARGET  COUNT,  NTARX , Y 

MAKE  NUMBER  OF  TARGET,  NCX,Y,  WITH  RXT , YT 
DETERMINE  TARGET  TYPE,  LX,Y 


COMPARE:  PKX , Y( LX , LY)  >  RX,Y,  TARGET  IS  HIT 
PKX , Y( LX , LY)  <  RX , Y ,  MISS 
PKHX , Y(LX , LY)  >  RX,Y,  HIT  BY  AIMED  FIRE 
PKHX , Y( LX , LY)  <  RX,Y,  HIT  BY  AIMED  FIRE 


1100  CONTINUE 

DO  1180  M=1 ,NGT 
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IF(RN(M).LE.(REAL(NTX)/NTT))  THEN 
X  SHOOTS 
NSX=NSX+1 

IF(NSX.GT.NTX)  GO  TO  1120 

IF(XSTAT(NSX) .EQ.O.O)  GO  TO  1180 

LX=NQX(NSX) 

IF(PKSX(LX) .LT.RXS(NSX) )  THEN 
NTARY=NTARY+1 
NCY=RYT ( NT ARY) *NTY+ 1 
LY=NQY (NCY) 

IF(PKX(LX,LY) .GE.RX(NSX) )  THEN 

KILL  TARGET  NUMBER  NCY 

IF(YSTAT(NCY) .GT.0.0)  THEN 
YSTAT(NCY)=0.0 
IF(PKHX(LX,LY) .GE.RX(NSX)) 
NDHY ( LY ) =NDHY ( LY ) + 1 
KHX (LX , LY) =KHX (LX , LY) +; 

ND AY ( LY ) =NDAY ( LY ) + 1 
KAX ( LX , LY ) =KAX ( LX , LY ) +  ] 
END  IF 
END  IF 
END  IF 

GO  TO  1180 
ELSE 

XSTAT (NSX) =0 . 0 
NDSX(LX)=NDSX(LX)+1 
GO  TO  1180 
END  IF 
ELSE 


Y  SHOOTS 


NSY=NSY+1 

IF(NSY.GT.NTY)  GO  TO 

IF(YSTAT(NSY) .EQ.O.O)  GO  TO 

LY=NQY(NSY) 

IF(PKSY(LY) . LT.RYS(NSY) )  THEN 
NTARX=NTARX+ 1 
NCX=RXT ( NT ARX ) *NTX+ 1 
LX=NQX(NCX) 

IF(PKY(LY,LX) .GE.RY(NSY)) 


1150 

1180 


THEN 


KILL  TARGET  NUMBER  NCX 


IF (XSTAT(NCX) .GT.0.0)  THEN 
XSTAT(NCX)=0 .0 
IF(PKHY(LY,LX) . GE .RY(NSY) ) 
NDHX ( LX ) =NDHX ( LX ) + 1 

ELSEKHY(LY-LX)=KHY(LY,LX)+1 

NDAX ( LX ) =NDAX ( LX ) + 1 
KAY ( LY , LX ) =KAY ( LY , LX ) + 1 
ENDIF 
ENDIF 
ENDIF 

GO  TO  1180 

ELSE 

YSTAT(NSY)=0.0 
NDS Y ( LY ) =NDSY ( LY ) + 1 
GO  TO  1180 
ENDIF 
ENDIF 


THEN 


THEN 


ALL  TYPE  X  COMBATANTS  HAVE  FINISHED  FIRING 
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CONTINUE 
KST0PX=1 
GO  TO  1170 


ALL  TYPE  Y  COMBATANTS  HAVE  FINISHED  FIRING 


CONTINUE 

KSTOPY=l 

EVERY  ONE  HAS  FINISHED  FIRING 
CONTINUE 

IF(KSTOPX,EQ . 1  .AND.  KST0PY.EQ.1) 
CONTINUE 

CALCULATE  THE  KILL  RATIOS 


GO  TO  1200 


CONTINUE 

DO  1290  J=1 ,NTYPX 
DO  1250  1=1 ,NTYPY 
NUMD=NX ( J ) *NY ( I ) 

IF  (NUMD.GE.l)  THEN 

DDD ( J , I ) =REAL ( KHX ( J , I ) *NTY ) / REAL (NUMD ) 
ELSE 

DDD(J,I)=0.0 
END  IF 

DDDAV ( J , I ) =DDDAV (  J ,  I ) +DDD ( J ,  I ) 

NUMB=NY ( I ) *NX ( J ) 

IF  (NUMB .GE . 1 )  THEN 

BBB(I , J)=REAL(KHY(I , J) *NTX) /REAL (NUMB) 
ELSE 

BBB(I , J)=0.0 
END  IF 

BBBAV ( I , J ) =BBBAV ( I ,  J ) +BBB ( I , J ) 

NUMC=NX ( J ) *NY ( I ) 

IF  (NUMC.GE.l)  THEN 

CCC( J, I )=REAL(KAX( J , I ) )/REAL(NUMC) 

ELSE 

CCC(J,I)=0.0 
END  IF 

CCCAV ( J ,  I ) =CCCAV ( J , I ) +CCC (  J ,  I ) 
NUMA=NY(l)*NX(J) 

IF  (NUMA.GE . 1 )  THEN 

AAA (I , J)=REAL(KAY(I , J) )/REAL(NUMA) 

ELSE 

AAA( I , J)=0 . 0 
END  IF 

AAAAV ( I , J ) =AAAAV ( I , J ) + AAA ( I ,  J  ) 

CONTINUE 

TOTAL  KILLED  BY  ALL  MECHANISMS 


NDX ( J ) =NDHX ( J ) +NDAX ( J ) +NDS  X ( J ) 
NDY ( JJ  =NDHY ( J ) +NDAY ( J ) +NDS Y ( J ) 
1290  CONTINUE 

IF ( IANS. EQ. I YES) THEN 
WRITE (6, 9560)  NSTEP , 

END  IF 


60)  NSTEP, ( (KHX (I ,J),J=1,3),I=1,3) 


FORMAT (/. 1  TIMESTEP' ,  14,/, 

X  10X/  AIMED  FIRE  KILLS  OF  Y"S  BY  X  FORCE 
X  .//V  9(1, 1)  D(l,2)  P ( 1 , 3 )  D(2,l)  D 


,//,'  D(1 , 1 )  D(1 , 2)  D( 1 ,3)  D(2 , 1 )  D(2,2)  D(2,3)', 

1  D(3,l)  D(3 , 2)  D(3 , 3) ' ,/ , IX, 

if(ians.eq.iyes)then 

WRITE(6,9665)  ( (KHY(I , J) , J=1 , 3) , 1=1 , 3) 


WRITE(6,9665)  ((KHY(I,J) ,J=1,3) ,1=1,3) 

END  IF 

9665  FORMAKIOX'  AIMED  FIRE  KILLS  OF  X"S  BY  Y  FORCE  1  , 

X  //,'  B(l.l)  B( 1 , 2)  B( 1 , 3)  B(2 , 1 )  B(2,2)  B(2,3) 
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Ln 


8 


[v 


9662 


9668 


9651 


C 

C 

C 


9650 


v  '  B(l4)/sB(3'2>  B(3,3)',/,1X, 

a  y lo , / ) 

IF ( IANS .EQ.IYES) THEN 

WRITER, 9682)  ((KAX(I,J)  ,J=1,3)  ,1=1,3) 

FORMAT (10X. '  AREA  FIRE  KILLS  OF  Y"S  BY  X  FORCE 

I  /A:  S&Ij  c<2'2>  c(2'i)' 

IF( IANS. EQ.IYES) THEN  ' 

WRITE(6,9§68)  ( (KAY (I , J) , J=1 , 3) , 1=1 ,3) 

FORMAT(10X  '  AREA  FIRE  KILLS  OF  X"S  BY  Y  FORCE 

l  /L'  iiy)  $;!)..?&»  M2'2> 

IF ( IANS .EO. IYES) THEN 
WRITE (6, 9700) 

So III > % • !I I ? (  NSTEP ' NTXS ' NTX ' NX  < NDX - NTYS , NTY , NY , NDY 

WRITE \1 ,9651 ) 

END  IF 

Fiw^s.Eo:!iES)i5LROWS  IH  DESCEND™° ■> 

88f{?:!SI8 

WRITE  7 ,9650  NSTEP,  KAX  I  J  J=1  3  1=1  3( 

2S^7'965°5  nstep , ( (kay 1 1 ; j) ; j=i ; 3 j ; i=i ; I > 

END  I  r 

F0RMAT(I4,9I8) 

CALCULATE  TOTAL  LOSSES  FOR  ALL  INTERACTIONS 


1500 


CONTINUE 
NDTX=0 
NDTY=0 

DO  1510  J=1 , 3 
NDTX=NDTX+NDX ( J ) 

NDTY=NDTY+NDY ( J ) 

1510  CONTINUE 

IF ( IANS. EQ. IYES )THEN 

WRITE (6,9670)  NDTX , NRRR , NDHX , NDAX , NDSX , 

XrmTTC.^  NDTY,NSSS,NDHY,NDAY,NDSY 

WRITE ( 7 , 9670 )  NDTX , NRRR , NDHX , NDAX , NDSX , 

„  _  -  X  ...  *  „  ,  NDTY ,  NSSS ,  NDHY ,  NDAY ,  NDSY 

9670  FORMAT (  _ TOTAL^  FORCE_  L0SS_ WAS1 ,  14 , 1 ,  RESUPPLY  WAS  s 1 , 314 , / , 


X 

X 

X 

X 

X 

X 

X 


AIMED  FIRE  LOSSES 


’,314,/, 
1 ,314,/, 
1 ,314,/, 


C 

C 

C 


ENDIF 


AREA  FIRE  LOSSES: 

SELF  LOSSES:  , 

TOTAL  Y  FORCE  LOSS  WAS1  ,14'’’  , 
AIMED  FIRE  LOSSES : 1 , Sl4,  / , 
AREA  FIRE  LOSSES 314, / , 
SELF  LOSSES: ’ ,314$ 


RESUPPLY  WAS: 1 ,314,/, 


ADD  IN  RESUPPLIED  UNITS 


1600  CONTINUE 

DO  1650  K=1 , 3 

IF(K.LE.NTYPX  .AND.  NRRR ( K ) . NE . 0 )  THEN 
NADDX=NRRR { K )  /oh 

NX ( K ) =NX ( K ) +NADDX 
DO  1620  1=1 ,NADDX 
NTX=NTX+ 1 
XPOS ( NTX ) =XLOC ( K ) 

NOX(NTX)=K 
XSTAT (NTX)=1 .0 
1620  CONTINUE 
ENDIF 

IF(K.LE.NTYPY  .AND.  NSSS(K) .NE.O)  THEN 
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NADDY=NSSS(K) 

NY ( K ) =NY ( K ) +NADDY 
DO  1640  1=1 ,NADDY 
NTY=NTY+1 
YPOS (NTY)=YL0C(K) 
NQY(NTY)=K 
YSTAT(NTY)=1 .0 
1640  CONTINUE 
END  IF 

1650  CONTINUE 


CLOSE  OUT  DESTROYED  UNITS 


1700 


1710 


1740 


1780 


CONTINUE 
DO  1710  J=1 , 3 
NX( J)=0 
NY(J)=0 
CONTINUE 
IX=0 

DO  1740  1=1, NTX 
IF(XSTAT(I) .NE.0.0)  THEN 
IX=IX+1 


NQX(IX)=NQX(I) 
XPOS(IX)=XPOS(I) 
XSTAT(IX)=XSTAT( 
K=NQX( I ) 
NX(K)=NX(K)+1 
END  IF 
CONTINUE 
NTX=IX 


I) 


IY=0 

DO  1780  1=1 ,NTY 
IF(YSTAT(I) .NE.0.0)  THEN 
IY=IY+1 


NOY ( IY ) =NQY ( I ) 
YPOS(IY)=YPOS(I) 

YSTAT?IY)=YSTAT(I) 

K=NQY ( I ) 

NY (K)=NY(K)+1 


END  IF 
CONTINUE 
NTY=IY 


PRINT  STATUS  FOR  EACH  TIME  STEP  IF  REQUESTED 
IF ( IANS . EO • IYES )THEN 

WRxTE(7 , 9/10 )  NSTEP , NTXS , NTX , NX , NDX , NTYS , NTY , NY , NDY 
WRITE (6, 9710)  NSTEP , NTXS , NTX , NX , NRRR , NTYS , NTY , NY , NSSS 
END  IF 


ACCUMILATE  STATISTICAL  VARIABLES  FOR  EACH  TIME  STEP 


I=NSTEP 


SNTX 

SNX1 

SNX2 

SNX3 

SNTY 

SNY1 

SNY2 

SNY3 

SSNTX 

SSNX1 

SSNX2 

SSNX3 

SSNTY 

SSNY1 


SSNY2 (I 
SSNY3 ( I 


SNTX( 
SNX1I 
SNX2  ( 
SNX3  ( 
SNTY  ( 
SNY1( 
SNY2( 
SNY3  ( 
=55^ 
=SSNX1 I 
=SSNX2  < 
=SSNX3 ( 
=SSNTY ( 
=SSNY1 ' 
=SSNY2 I 
=SSNY3i 


'I ' 

'V 

^x'( 


+NTX 
+NX  ( 1 ) 

+NX ( 2 ) 

+NX ( 3 ) 

+NTY 

+NY( 1 ) 

+NY(2) 

+NY(3) 

I)+NTX**2 

I)+NX(1)**2 

I)+NX(2)**2 

I)+NX(3  **2 

I)+NTY**2 

I)+NY(1)**2 

I)+NY(2)**2 

I)+NY(3)’I',I'2 
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TERMINATE  IF  FORCES  ARE  ATTRITTED  BELOW  SPECIFIED  LEVEL 
I F ( NTX . LE . NXSTOP  .OR.  NTY.LE.NYSTOP)  GO  TO  9120 
TERMINATE  ON  NUMBER  OF  CYCLES 


IF(NSTEP.LT.NSTOP) 
9120  CONTINUE 

STEPS=REAL(NSTEP) 


GO  TO  1000 


CALCULATE  THE  AVERAGE  ATTRITION  COEFFICIENTS  PER  REPLICATION 


DO  9100  J=1 ,3 
DO  9100  1=1,3 

AAAAV ( I , J ) =AAAAV ( I , J ) /STEPS 
BBBAV ( I ,  J  J =BBBAV ( I ,  J ) /STEPS 
CCCAV ( I ,  J  =CCCAV  I , J ) / STEPS 
DDDAV(I , J)=DDDAV(I , J) /STEPS 
S AAAAV ( I , J ) =SAAAAV ( I ,  J ) +AAAAV ( I , J 1 
SBBBAV ( I , J > =SBBBAV (I , J) +BBBAV ( I , J 
SCCCAV (I , J) =SCCCAV ( I ,  J) +CCCAV { I ,  J 
_ , SDDDAV (I , J) =SDDDAV (I , J) +DDDAV ( I , J 
9100  CONTINUE  v  '  1 

CSTEP (N ) =NSTEP 


ACCUMILATE  THE  STATISTICAL  VARIABLES  FOR  THE  NUMBER  OF  TIME  STEPS 

SSTEP=SSTEP+CSTEP (N ) 

SSSTEP=SSSTEP+(CSTEP(N) **2) 

WRITE  THE  NUMBER  OF  REPLICATIONS  AND  THE  NUMBER  OF  STEPS  TO  FILE 


WRITE (12 , 8510 )N, CSTEP (N) 
8510  FORMAT ( 13 , F4 . 0 )  ' 


KEEP  TRACK  OF  MINIMUM  NUMBER  OF  STEPS  FOR  ALL  REPLICATIONS 


IF (MINST .GT . CSTEP (N) ) THEN 
MINST=CSTEP (N) 

ENDiF 

4000  CONTINUE 


CALCULATE  AVERAGE/STANDARD  DEV. OF  FORCE  LEVELS  FOR  EACH  TIME  STEP 


B=REAL(L) 

RD=1 . 0/B 

RDD=1.0/(B*(B-1)) 
DO  4400  1=1, MINST 


ANTX 

ANX1 

ANX2 

ANX3 

ANTY 

ANY1 


!=SNTX( 
=SNX1 
=SNX2 
=SNX3 
=SNTY 
=SNY1 ( 


ANY2(I)=SNY2I 


ANY3u 

SDNTXI 

SDNX1 < 

SDNX2( 

SDNX3 

SDNTY 

SDNY1 

SDNY2 

SDNY3 


>=SNY3(I)*RI 
I)=SQRT(ABSI 
I  )=SQRT(ABS( 
I  )=S$RT (ABS ( 
[  )=SQRT (ABS ( 
[  )=SQRT (ABS ( 
C )=SQRT( ABS ( 
[ ) =SQRT (ABS ( 
[  )=SoRT( ABS ( 


iB*SSNTX(i; 
B*SSNX1 ( I 
B*SSNX2 ( I 
B*SSNX3 ( I 
B*SSNTY(I 
B*SSNY1  I 
B*SSNY2 ( I 
B*SSNY3 ( I ) 


■SNTX( I 
-SNX1 ( I 
-SNX2 ( I 
-SNX3 ( I 
-SNTY ( I 
-SNY1 ( I ' 
-SNY2 ( I ' 
■SNY3 { I ' 


i**2)*RDD; 
l**2)*RDD] 
>**2)*RDD 
**2 >*RDD) 
**2>*RDD) 
**2)*RDD 
**2)*RDD) 
**2 ) *RDD ) 


?So5U^IE«Igg  Ki  fgg  HINUS  0NE  S7mmm  AVIATION 


•  */  V  V  V  /  /  v  V  /  /  ’ 
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PSNTX ( I ) =ANTX ( I ) +SDNTX ( I ) 

PSNX1 ? I ) =ANX1 ? I ) +SDNX1 ? I ) 

PSNX2 ( I ) =ANX2 ( I ) +SDNX2 ? I 
PSNX3 ( I )=ANX3 ( I ) +SDNX3 ? I ) 

PSNTY ( I ) =ANTY ( I ) +SDNTY ( I ) 

PSNY1 ( I )=ANY1 ( I ) +SDMY1 ( I 
PSNY2 ( I ) =ANY2 ( I ) +SDNY2 ( I ) 

PSNY3 ( I ) =ANY3 ( I ) +SDNY3 ( I ) 

TSNTX? I )=ANTX( I ) - (SDNTX( I ) ) 

TSNX1 ( I ) =ANX1 ( I ) -SDNX1 ( I ) 

TSNX2 ( I ) =ANX2 ( I ) -SDNX2 ( I ) 

TSNX3 ( I )=ANX3 ( I ) -SDNX3 ( I ) 

TSNTY ( I ) =ANTY ( I ) - SDNTY ( I ) 

TSNY1 ( I )=ANY1 ( I ) -SDNY1 ( I ) 

TSNY2  I ) =ANY2 ( I ) - SDNY2  I 
TSNY3 ( I ) =ANY3 ( I ) -SDNY3 ( I ) 

CALCULATE  THE  AVERAGE  ATTRITION  COEFFICIENTS  OVER  ALL  REPLICATIONS 


DO  9200  J=1 , 3 


DO  9200  K=1 , 3 

AAAAV ( K , J ) =SAAAAV ( K , J ) 
BBBAV ( K , J ) =SBBBAV ( K , J ) 
CCCAV(K, J  =SCCCAV(K, J 
DDDAV ( K , J ) =SDDDAV ( K , J ) 
CONTINUE 


*RD 

*RD 

*RD 

*RD 


PRINT  STATISTICS  TO  FILES  8,9,10,11  FOR  PLOTTING 


WRITE(8,8100)  I , ANTX ( I ) , ANTY ( I ) , PSNTX ( I ) , PSNTY ( I ) , TSNTX ( I ) , 
XTSNTY (I ) 

WRITE (9, 8100)  I , ANX1 ( I ) , ANY1 ( I ) , PSNX1 (I ) , PSNY1 ( I ) , TSNX1 ( I ) , 
XTSNYl(I) 

WRITE { 10 , 8100 )  I , ANX2 ( I ) , ANY2 ( I ) , PSNX2 ( I ) , PSNY2 ( I ) ,TSNX2(I) 
XTSNY2 (i ) 

WRITE? 11 ,8100)  I , ANX3 ( I ) , ANY3 ( I ) , PSNX3 ( I ) , PSNY3 ( I ) , TSNX3 ( I ) 
XTSNY3 (I ) 

8100  FORMAT? 13 , 6 (F8. 2 , IX) ) 


WRITE  AVG'S  AND  STANDARD  DEVIATIONS  TO  FILE  FOR  TABLES 


WRITE (13,8150)1 , ANTX (I ) , SDNTX( I ) , ANTY(I ) , SDNTY ( I 
WRITE (14,8150)1 , ANX1 ( I ) , SDNX1 ( I ) , ANY1 ( I ) , SDNY1 ( I 
WRITE (15,8150)1 , ANX2 ( I ) , SDNX2 ( I ) , ANY2 ( I ) ,SDNY2(I 
WRITE (16,8150)1 , ANX3 ( I ) , SDNX3 ( I ) , ANY3 ( I ) , SDNY3 ( I 
8150  FORMAT ( I 3, 4 (F8. 2, IX) ) 

4400  CONTINUE 


CALCULTE  THE  AVG  AND  STANDARD  DEVIATION  FOR  THE  #  OF  NSTEPS  PER  RUN 


AVSTEP=SSTEP*RD 

SDSTEP=SQRT (ABS ( (B*SSSTEP-SSTEP**2 )*RDD) ) 
WRITE (12,8500 )MINST , AVSTEP , SDSTEP 
8500  FORMAT? 13 ,2F7 .2) 


EXPLAIN  CAUSE  OF  TERMINATION 


9000  CONTINUE 


WRITE ( 
9656  FORMAT? 
WRITE? 
WRITE? 
9655  FORMAT ( 
WRITE? 
WRITE? 
9690  FORMAT? 
X  10X, 1 
X 
X 


7,9656) 

'  a  b  r  d  BY  ROWS  IN  DESCENDING  ORDER  ' ) 

?  AAAAV  ?I,J),J=1,3),I  =  1,3) 

?BBBAV?I,J),J=1,3),I=1,3) 

8.5) 

6,9695)  ? (AAAAV? I, J) ,J=1,3) ,1=1,3) 

6,9697)  (?BBBAV?I,J) ,J=1,3) ,1=1,3) 

/, '  TIMESTEP1 ,14,/, 

AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE  (  X  E+03)',//, 
'  CC?1,1)  CC ( 1 , 2 )  CC ? 1 , 3 )  CC(2 , 1 )  CC(2,2)  CC( 2,3)', 

•  CC?3 , 1 )  CC  v  2)  CC(3 , 3) ' ,/, 


-w 

7,9655 
7,9655)  ( 
1  AV  1 ,9F 


nrm 


X  3P ,9F8 .  5 ,  / ) 

WRITE(7,9655)  ( (CCCAV(I , J) , J=1 , 3) , 1=1 , 3) 

WRITE (7, 9655)  ( (DDDAV( I , J) , J=1 , 3) , 1=1 , 3 
WRITE  6,9698)  CCCAV  I , J) , J=1 , 3  , 1=1 , 3 
WRITE  6,9699)  ( (DDDAV(I , J) , J=1 , 3) , 1=1 , 3 ) 

9695  FORMAT (10X, 

X  1  AVERAGE  ATTRITION  COEFFICIENTS  FOR  Y  FORCE 
X  •  AA(1 , 1 )  AA( 1,2)  AA(  1 , 3)  AA(2,1)  AA(2,2)  AA(2,3) 

X  *  AA(3 , 1 )  AA(3 , 2)  AA(3,3)' ,/,lX, 

X  9F8 . 5 , / ) 

9697  FORMAT (10X, 

X  1  AVERAGE  ATTRITION  COEFFICIENTS  FOR  Y  FORCE 
X  '  BB(1 , 1 )  BB(1 , 2)  BB ( 1 , 3 )  BB(2,1)  BB(2,2)  BB(2,3) 

X  '  BB(3,1)  BB(3,2)  BB(3 , 3) 1 , / , IX, 

X  9F8 . 5 , / ) 

9698  FORMAT(10X, 

X  '  AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE 
X  1  CC(1,1)  CC{1 ,2)  CC( 1 , 3)  CC(2 , 1 )  CC(2,2)  CC(2,3)' 

X  1  CC(3,1  CC(3,2)  CC(3 , 3 ) 1 , / , IX, 

X  9F8.5,/) 

9699  FORMAT (10X, 

X  1  AVERAGE  ATTRITION  COEFFICIENTS  FOR  X  FORCE 
X  1  DD( 1 , 1 )  DD( 1 , 2)  DD( 1 , 3)  DD(2,1)  DD(2,2)  DD(2,3)' 

X  1  DD(3 , 1 j  DD(3 , 2 )  DD  3 , 3 ) 1 , / , IX , 

X  9F8. 5 , / ) 


THEN 

NSTEP , NTY , NTYS 
NSTEP,NTY,NTYS 


I F ( NTX . LT . NXSTOP )  THEN 

WRITE  (6,9620)  NSTEP , NTX, NTXS 
WRITE  (7 ,9620)  NSTEP, NTX, NTXS 
ENDIF 

IF (NTY.LT.NYSTOP)  THEN 

WRITE  (6,9630)  NSTEP, NTY, NTYS 
WRITE  (7, 9630)  NSTEP, NTY, NTYS 
ENDIF 

IF  (NSTEP  .EQ.  NSTOP)  THEN 
WRITE  (6,9640)  NSTOP 
WRITE  (7,9640)  NSTOP 
ENDIF 

9620  FORMAT ( 1  STOPPED  IN  STEP  ',15,'  BY  FORCE  X',/, 

X  '  CURRENT  X  STRENGTH  =',I5,',  INITIAL  X  STRENGTH1 

9630  FORMAT ('  STOPPED  IN  STEP  ',15,'  BY  FORCE  Y',/, 

X  1  CURRENT  Y  STRENGTH  = ' , 15 , ' ,  INITIAL  Y  STRENGTH' 

9640  FORMAT ( 1  STOPPED  BY  PROGRAM  AT  NSTOP  =  ’,15,//) 

END 


,15,//) 

,15,//) 


i 


APPENDIX  C 

EXEC  FILE  FOR  M/VV  PROGRAM 


&TRACE  ON 

GLOBAL  TXTLIB  VFORTLIB  NONIMSL  IMSLSP  CMSLIB 

GLOBAL  LOADLIB  VFLODLIB 

FILEDEF  4  TERM 

FILEDEF  5  DISK  INMWM  DATA  A 

FILEDEF  6  DISK  MWM  STLIST  A 

FILEDEF  8  DISK  TOTAL  DATA  A 

FILEDEF  9  DISK  FORI  DATA  A 

FILEDEF  10  DISK  FOR2  DATA  A 

FILEDEF  11  DISK  FOR3  DATA  A 

LOAD  MWM 

START  MAIN 

CP  MSG  0731P  JOB  MWM  COMPLETED 
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APPENDIX  D 

M/W  COMPUTER  PROGRAM 

PROGRAM  WRITTEN  TO  RUN  THE  GENERALIZED  LANCHESTER 
EQUATION  MODEL  DESIGNED  BY 
PAUL  H.  MOOSE  AND  JACK  M.  WOZENCRAFT 
1/08/87 

WRITTEN  TO  TEST  THE  MOOSE/WOZENCRAFT  THEORY 


t 


NTYPX 

NTYPY 

TNX(I) 

TNY(J) 

STNX,Y 

TTNX ,  Y 

AAA(J,I] 
BBB(  J , I 
CCC(I,  J 
DDD(I , J) 

SAAA(I) 

SBBB(I) 

SCCC(J) 

SDDD(J) 


ARLX(I) 

ARLY(J) 

SLX(I) 
SLY( J) 

TNDX(I) 

TNDY(J) 

NSTOP 

XSTOP 

YSTOP 


NUMBER  OF  TYPES  OF  X  COMBATANTS  <INPUT> 

NUMBER  OF  TYPES  OF  Y  COMBATANTS  <INPUT> 

NUMBER  OF  EACH  TYPE  OF  X  COMBATANTS 
NUMBER  OF  EACH  TYPE  OF  Y  COMBATANTS 
VARIBLE  TO  SUM  TOTAL  FORCE  LEVEL  PER  TIMESTEP 
TOTAL  FORCE  LEVEL  PER  TIMESTEP 

X(I)  AREA  FIRE  LOSSES  DUE  TO  Y(J) 
ATTRITION:  X(I)  AIMED  FIRE  LOSSES  DUE  TO  Y  J) 

ATTRITION:  Y(J)  ARERA  FIRE  LOSSES  DUE  TO  X(I) 

ATTRITION:  Y( J )  AIMED  FIRE  LOSSES  DUE  TO  X(I) 

ATTRITION:  X(I)  SUM  OF  THE  AREA  LOSSES 

SUM  OF  THE  AIMED  FIRE  LOSSES 
ATTRITION:  Y(J)  SUM  OF  THE  AREA  FIRE  LOSSES 
ATTRITION:  Y(J)  SUM  OF  THE  AIMED  FIRE  LOSSES 

X(I)  RESUPPLIED/TIMESTEP 
NUMBER  COMBATANT  Y(J)  RESUPPLIED/TIMESTEP 

SELF-ATTRITION  COEFFICIENT  FOR  X  FORCE 
SELF-ATTRITION  COEFFICIENT  FOR  Y  FORCE 

X(I)  LOSSES  TO  AIMED  FIRE  THIS  TIMESTEP 
Y(J)  LOSSES  TO  AIMED  FIRE  THIS  TIMESTEP 

X(I}  LOSSES  TO  AREA  FIRE  THIS  TIMESTEP 
Y(J)  LOSSES  TO  AREA  FIRE  THIS  TIMESTEP 


t(I)  SELF  LOSSES  THIS  TIMESTEP 
r(J)  SELF  LOSSES  THIS  TIMESTEP 


TOTAL  CHANGE  IN  X(I) 
TOTAL  CHANGE  IN  Y(J) 


FORCE  LEVEL  THIS  TIME  STEP 
FORCE  LEVEL  THIS  TIME  STEP 


MAXIMUM  NUMBER  OF  TIMESTEPS  <INPUT> 

X  FORCE  LEVEL  DEFINED  AS  THE  BREAK  PT. 
Y  FORCE  LEVEL  DEFINED  AS  THE  BREAK  PT, 


IMPLICIT  REAL*8  (D) 

CHARACTER*40  I HEAD 
CHARACTER* 1  IANS,IYES 

DIMENSION  TNX ( 3 ) , TNY ( 3 ) 
DIMENSION  NRRR(3) ,NSSS(3' 
DIMENSION  AMLX(3) ,AMLY(3! 
DIMENSION  ARLX(3) ,ARLY(3 


W'Vv 


DIMENSION  SLX(  3 ) , SLY ( 3 ) 
DIMENSION  TNDX(3) ,TNDY(3) 
DIMENSION  AAA (3, 3  ,SAAA(3) 
DIMENSION  BBB (3,3  ,SBBB(3 
DIMENSION  CCC(3,3),SCCC(3) 
DIMENSION  DDD(  3 ,  3 ) , SDDD( 3 ) 
DIMENSION  UUU ( 3 ) ,  VW ( 3 ) 

C 

C 

DATA  IYES/'Y'/ 


CALL  EXCMS ( 1 CLRSCRN ' ) 

READ  (5,9400)  I HEAD 
WRITE(4,9410) 

9410  FORMAT  '  DO  YOU  WANT  TO  CHANGE  THIS  HEADING?  (Y/N)  (A40)‘) 
WRITE (4, 9400)  I HEAD 
READ  (4,9405)  IANS 
9405  FORMAT(Al) 

IF( IANS . EQ . IYES )  READ  (4,9400)  I HEAD 

9400  FORMAT (A40) 

9401  FORMAT(Al) 

WRITE (6, 9400)  I HEAD 

WRITE (8, 9400)  I HE AD 

WRITE (9, 9400  I HEAD 

WRITE (10 ,9400)  I HEAD 

WRITE (11, 9400)  I HE AD 

READ  (5,9500)  NSTOP 

READ  (5,9500)  NTYPX , TNX , UUU , NRRR 

WRITE(6  9502) 

WRITE(6  9501 )  NTYPX , TNX . UUU , NRRR 
READ  5,9505)  ( (DDD(I , J) , J=1 , 3 ) , 1=1 , 3) 

WRITE  6,9521) 

WRITE (6 , 9520 )  ((DDD(I, J) , J=1 ,3) ,1=1,3) 

READ  (5,9505)  ( (CCC( I , J ) , J=i , 3 ) , 1=1 , 3 ) 

WRITE(6 ,9523) 

WRITE(6 , 9520 )  ( (CCC( I , J) , J=1 , 3 ) , 1=1 , 3 ) 

READ  5,9500  NTYPY,TNY,VW,NSSS 
WRITE (6  9503) 

WRITE (6 | 9501 )  NTYPY,TNY, WV,NSSS 
READ  (5,9505  ( (BBB( I , J) , J=1 , 3)  ,  1  =  1 , 3  ) 

WRITE (6  9522) 

WRITE (6; 9520  ((BBB(I, J) , J=l,3) ,1=1,3) 

READ  (5,9505)  ( (AAA( I ,J),J=1,3),I=1,3) 

WRITE (6  9524) 

WRITE (6 ' 9520)  ((AAA(I.J),J=1 ,3) ,1=1,3) 

9500  FORMAT (I5,3F5.0,3F7.3,3I5) 

9501  FORMAT (15, 3F7 . 2 , 3F7 . 3 , 3 1 5 ) 

9502  FORMAT (  1  #  X  TYPES  NUMBER  EACH  SELF-LOSS 

X  ’  RE -SUPPLY  1  ) 

9503  FORMAT (  1  #  Y  TYPES  NUMBER  EACH  SELF-LOSS 

X  1  RE-SUPPLY  ' ) 

9505  FORMAT (3F8. 5) 

9520  FORMAT (9F8. 5,/) 

9521  FORMAT ( / , '  KILL  PROBABILITIES  FOR  FORCE  X(I)  SHOOTS  Y(J)1,//, 

X  '  DD( 1 , 1 )  DD( 1 , 2)  DD  1,3  DD(2,1)  DD(2,2)  DD ( 2 , 3 ) ' , 

X  '  DD(3 , 1 )  DD[3,2)  DD(3,3) 1 ) 

9522  FORMAT (/ , 1  KILL  PROBABILITIES  FOR  FORCE  Y(I)  SHOOTS  X(J)',//, 

X  '  BB ( 1 , 1 )  BB( 1 , 2 )  BB  1,3  BB(2,1)  BB(2,2)  BB(2,3)', 

X  1  BB  3,1)  BB(3 , 2)  BB(3,3) ' ) 

9524  FORMAT (/ , 1  AREA  FIRE  COEFFICIENTS  FORCE  Y(I)  SHOOTS  X(J)  ',//, 
X  1  AA( 1 , 1 )  AA ( 1 , 2 )  AA ( 1 , 3 )  AA(2,1)  AA(2,2)  AA(i,3)', 

X  1  AA( 3  l)  AA( 3  2)  AA(3  3)') 

9523  FORMAT (/ , 1  AREA  ^IRE  COEFFICIENTS  FORCE  X(I)  SHOOTS  Y(J)  '  .11, 

X  1  CC(1,1)  CC( 1 , 2 )  CC ( 1 , 3 )  CC(2 , 1)  CC(2,2)  CC(2,3)‘, 

X  1  CC( 3 , 1 )  CC(3 , 2 )  CC ( 3 , 3  )  1  ) 


C 

C 

C 


INITIALIZE  THE  INTIAL  FORCES  LEVELS 


INTX=Q 

INTY=Q 

DO  1100  1*1 , 3 

INTX=INTX+TNX( I ) 

INTY=INTY+TNY(I ) 

1100  CONTINUE 

TTNX=INTX 
TTNY=INTY 
XSTOP=50 .0 
YSTOP=50 . 0 
MINST=99 

;  SETUP  A  FILE  FOR  TRACING  EACH  STEPS  RESULTS 

WRITE (6,9711) 

9711  FORMAT { 1  #  INTX  TTNX  TNX1  TNX2  TNX3  RX1  RX2  RX3 ' 

X  '  INTY  TTNY  TNY1  TNY2  TNY3  RY1  RY2  RY3 ' ) 

9712  FSSiSIl  I4,72(?l  INTY. TINY. TOT. NS5S 

WRITE (6 , 9700 )  * 

9700  FORMAT ( '  #  INTX  TTNX  TNX1  TNX2  TNX3 ' 

X  1  INTY  TTNY  TNY1  TNY2  TNYi 1 ) 

9710  PORSil452(lUf7E2iI"r'T™X'TN,!'INTY'TTNY’m 
:  END  IF 


THIS  IS  THE  BEGINNING  OF  THE  COMBAT  CYCLE 


DO  4000  N=1 , NSTOP 
NSTEP=N 


C  ZERO  STORAGE  FOR  THIS  TIMESTEP 

STNX=0 . 0 
STNY=0 . 0 

DO  1015  1=1 , NTYPX 
SAAA( I ) =0 . 0 
SBBB ( I ) =0 . 0 
1015  CONTINUE 

DO  1020  J  = 1 , NTYPY 
SCCC( J)=0.0 
SDDD( J)=0 . 0 
^1020  CONTINUE 

;  CALCULATE  THE  LOSSES  DUE  TO  AREA  FIRE  FOR  THIS  TIMESTEP 

DO  1110  1=1, NTYPX 
DO  1105  J=l, NTYPY 

J.!)*TNY(J)) 

1105  cONT¥NCuVjj=SCCClji  +  lcCc!l'ji*TNxilii 
1110  CONTINUE 

DO  1115  1=1, NTYPX 
,  ,  ARLX( I )*(TNX( I ) *SAAA ( I ) ) 

1115  CONTINUE 

DO  1120  J=l, NTYPY 

ARLY(J)=TNY(J)*SCCC( J) 

1120  CONTINUE 

CALCULATE  THE  LOSSES  DUE  TO  AIMED  FIRE  PER  TIMESTEP 

DO  1130  1=1, NTYPX 
DO  1125  J=l, NTYPY 

SBBB (I)=SB8B(I)+(BBB(J,I) *TNY ( J ) ) 

1125  CONT®  J^SDDD'J^'DDD'I’J^TNX'1^ 

1130  CONTINUE 

DO  1135  1=1, NTYPX 


non  non  nnnnnnn  nnnn  non  non  ono  noo 


AMLX( I )=(TNX( I )/TTNX) *SBBB( I ) 

1135  CONTINUE 

DO  1140  J*1 ,NTYPY 

AMLY ( J ) * ( TNY ( J ) /TTNY ) *SDDD ( J ) 

1140  CONTINUE 

CALCULATE  THE  SELF  LOSSES  PER  TIMESTEP 

DO  1142  1=1 ,NTYPX 

SLX ( I ) =TNX ( I ) *UUU ( I ) 

1142  CONTINUE 

DO  1145  J=1 , NTYPY 

SLY  ( J )  =TNY  ( J )  *VW  ( J ) 

1145  CONTINUE 

CALCULATE  THE  TOTAL  KILLED  BY  ALL  MECHANISMS 
DO  1150  1=1 , NTYPX 

TNDX( I )=-SLX( I ) -ARLX( I ) -AMLX( I )+NRRR( I ) 
1150  CONTINUE 

DO  1155  J=l, NTYPY 

TNDY ( J )=-SLY ( J ) -ARLY ( J ) - AMLY( J ) +NSSS ( J ) 
1155  CONTINUE 

CALCULATE  THE  FORCE  LEVELS  AFTER  THIS  TIMESTEP 

DO  1160  1=1, NTYPX 

TNX( I )=TNX( I ) +TNDX( I ) 

1160  CONTINUE 

DO  1165  J=l, NTYPY 

TNY ( J ) =TNY ( J ) +TNDY ( J ) 

1165  CONTINUE 

CALCUTE  THE  TOTAL  X,Y  FORCE  LEVELS 

DO  1170  1=1, NTYPX 
STNX=STNX+TNX ( I ) 

1170  CONTINUE 

TTNX=STNX 
DO  1175  J=l, NTYPY 
STNY=STNY+TNY ( J ) 

1175  CONTINUE 

TTNY=STNY 


PRINT  STATUS  FOR  EACH  TIME  STEP 

WRITE ( 6 , 67 10 )  NSTEP. INTX , TTNX , TNX , INTY , TTNY , TNY 
6710  FORMAT (14,2(15, 4F7 .  2 ) ,  / ) 


PRINT  VALUES  TO  FILES  8,9,10,11  FOR  PLOTTING 


WRIT1(8 , 8100 )  N, TTNX, TTNY 
WRITE(9 ,8100 )  N , TNX ( 1 ) , TNY ( 1 ) 

WRITE(i6,8106)  N , TNX ( 2 ) , TNY ( 2 ) 

WRITE(ll ,8100)  N,TNX(3) ,TNY(3) 

8100  FORMAT ( I 3, 2 (F8. 2, IX) ) 

TERMINATE  IF  FORCES  ARE  ATTRITTED  BELOW  SPECIFIED  LEVEL 
I F ( TTNX . LE . XSTOP  . OR . TTNY . LE . YSTOP )  GO  TO  5000 
RETURN  AND  CALULATE  FORCE  LEVELS  FOR  ANOTHER  TIMESTEP 
4000  CONTINUE 


'4 


noon 


EXPLAIN  CAUSE  OF  TERMINATION 


5000  CONTINUE 


IF(TTNX.LT.XSTOP)  THEN 

WRITE  (6,9620)  NSTEP , TTNX , INTX 
END  IF 

IF(TTNY.LT.YSTOP)  THEN 

WRITE  (6,9630)  NSTEP , TTNY , INTY 
END  IF 

IF  (NSTEP  ,EQ.  NSTOP )  THEN 
WRITE  (6,9640)  NSTOP 
Ef'DIF 

9620  FORMAT ( '  STOPPED  IN  STEP  ',15,'  BY  FORCE  X',/, 

X  1  CURRENT  X  STRENGTH  =',F5.2,',  INITIAL  X  STRENGTH 1  , 1 5 ,// ) 

9630  FORMAT ( '  STOPPED  IN  STEP  ',15,'  BY  FORCE  Y',/, 

X  1  CURRENT  Y  STRENGTH  =',F5.2,',  INITIAL  Y  STRENGTH ', 1 5 ,// ) 

9640  FORMAT ('  STOPPED  BY  PROGRAM  AT  NSTOP  =  ',15,//) 

END 


APPENDIX  E 

FORCE  LEVEL  TRAJECTORIES 


This  Appendix  is  designed  to  provide  the  reader  with  a  complete  set  ol  all  of  the 
results  of  this  experiment  T  he  input  data  sets  lor  each  model  are  presented  in  tables 
and  then  each  of  the  force  lex  el  trajectories  is  shown  graphicallx 


TABLE  1 

SAMPLE  INPUT  DATA  SET  TOR  C  S  MODEL 


<>l  16  S7  TEST  OF  (2X2).  1  UNSTABLE  ROOT  S 
NSTEP 

NTVPX  NTX1  N  1X2  NTX3  loci  loc2  loc3  u 


11 


d? 


2! 


Ml 

"21 

'1 


"}  1 2  *4  i  3 

Jr! 

"12  "13 

C  •»  i  C  1  3 

"32  "33 


1  u2  u3  rl  r2  r  3 


NTYPY  NTYl 
b.  bn 
bU  bU 

b?S  b^-s 
a 


b“ 

Ml 

til 


NTY2  NTY3  loci  loc2  loc3 


V1  v2 


3  S1 


s*>  Si 


'Vi 

.  1 2  d  P 

1  Jis  '*">3 

1  a32  aO 
DSEED 


DSL  ED 


DSEED 


I  ABLE  2 

SAMPLE  INPUT  DATA  SET  TOR  M  W  MODE 


"1  16  S’  I  ES  I  OF  i  2X2).  1  L  N  STAB  EL  ROOT  S 
\SI  EP 

NIYPX  N  1X1  VEX 2  NTX3  u ,  u,  u3  r,  r  3 

'".Ui’Jii 

"1 1  "1 2  "1 3 

c  ■*  j  C-O 

"M  "3  2  "33 

N1YPY  NTYl  NTY2  NTY3 

hi  h:  hi 


v|  v2  v3  "t  s2  s3 


TABLE  3 

INPUT  DATA  SET  FOR  C  S  MODEL  CASE  I 


01  05  87  TEST  OF  (1X1).  NO  ROOTS 


99 

1  600  0  0  1.0  1.0  1.0 

0.06000  0.<XXXX)  0.00000 
0.0OO00  0.00000  0. 00000 
((.(MMX XJ  O.(XXMH)  OtXHMM) 

0.00020  0.00000  OtXMMX) 

0.<XXXX)  O.OiXXX)  0.00000 
0.00000  0.00000  O.(XKKX) 

1  400  0  o  1.0  1.0  1.0 

0.0 1000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0  tXX>20  O.(KHKX)  O.WXXX) 

OAMXXX)  O.(XHMK)  O/XXXK) 

O.OO(XX)  O.(XXMX)  04XXMX) 

3452 15789.  DO  45635761.  DO 


0.040  O.(XX)  0.000  28  0  0 


0.010  0.000  O.(XX)  16  0  0 


89342761.  DO 


TABLE  4 

INPUT  DATA  SET  FOR  M  W  MODEL  CASE  1 


0105  87  TEST  OF  (1X1).  NO  ROOTS 
99 


1  600  0  0  1.0  1.0 

0.06000  0.00OO0  0.00000 
0.00000  0.00000  0.00000 
o .<xxxx)  o.ooooo  0.00000 

<MXX>20  0.00<XX)  0.(XXXH) 
0.00000  0.00000  O.(XXXX) 
0.00000  0.00000  O.(XXKX) 

1  400  o  o  1.0  1.0 

0.01000  0.00000  O.(XXXX) 
O.(XHXX)  0.0(HXX)  O.(XXXK) 
(MXXXX)  0.00000  O.OOCMXj 
0.<XX)20  0.00000  0.00000 
0.00000  0.00000  O.iXXXX) 
0.00000  0.(XXXX)  O.(XXXX) 


1.0  0.040  0.000  0.000  28 


1.0  0.010  0.000  0.000  16 


0 


0 


0 


0 


TOTAL  FORCE  LEVEL  TRAJECTORY 


TOTAL  X  FORCE  LEVEL 


CASE  1  (1X1)  NO  ROOTS 


Figure  E.l  Total  Force  Level  Trajectory  For  Case  One. 


X  FORCE  LEVEL 

200  250  300  350  400  450  500  550  600  650  700  750  BOO 


X  FORCE  LEVEL  TRAJECTORY 


TIME 


CASE  1  (1X1)  NO  ROOTS 


Figure  E.2  X  Force  Level  Trajectory  Over  Time  For  Case  One. 
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TABLE  5 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  2 


01  05/87  TEST  OF  (1X1), 
99 

1  600  0  0  1.0  1.0 
0.01000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
1  500  0  0  1.0  1.0 

0.01000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
345215789. DO 


1  STABLE  ROOT  AT  200.200S 
1.0  0.040  0.000  0.000  18  0 


1.0  0.060  0.000  0.000  22  0 


89342761. DO 


45635761. DO 


0 


0 


TABLE  6 

INPUT  DATA  SET  FOR  M  W  MODEL  CASE  2 


01  05/87  TEST  OF  (1X1),  1  STABLE 

1  600.  0.  0.  0.040  0.000  0.000 
0.00877  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00018  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  O.00000  0.00000 

1  500.  0.  0.  0.060  0.000  0.000 
0.00874  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00017  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 


ROOT  AT  200.200S 
18  0  0 


22  0  0 


TOTAL  FORCE  LEVEL  TRAJECTORY 


50  100  150  200  250  300  350  400  450  500  560  600  650  TOO 


TOTAL  X  FORCE  LEVEL 

CASE  2  (1X1)  1  STABLE  ROOT 

Figure  E.4  Total  Force  Level  Trajectory  For  Case  Two 
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X  FORCE  LEVEL  TRAJECTORY 


Y  FORCE  LEVEL  TRAJECTORY 

o 


TIME 

CASE  2  (1X1)  1  STABLE  ROOT 


Figure  E.6  Y  Force  Level  Trajectory  Over  Time  For  Case  Two. 
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TABLE  7 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  3 


01  06/87  TEST  OF  (1X1),  1  UNSTABLE  ROOT  AT  200.200S 
99 

1  500  0  0  1.0  1.0  1.0  0.010  0.000  0.000  22  0  0 

0.06000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.000 20  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

1  600  0  0  1.0  1.0  1.0  0.010  0.000  0.000  18  0  0 

0.04000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

345215789. DO  45635761. DO  89342761. DO 


TABLE  8 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  3 


01  06'87  TEST  OF  (1X1),  1  UNSTABLE  ROOT  AT  200.200S 
99 

1  500.  0.  0.  0.010  0.000  0.000  22  0  0 

0.05129  0.00000  0.00000 
O.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00018  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
1  600.  0.  0.  0.010  0.000  0.000  18  0  0 

0.03460  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00017  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 


TOTAL  Y  FORCE  LEVEL 


TOTAL  FORCE  LEVEL  TRAJECTORY 

o 


SO  100  ISO  200  250  300  360  400  460  600  650  600  660  700 

TOTAL  X  FORCE  LEVEL 

CASE  3  (1X1)  1  UNSTABLE  ROOT 


Figure  E.7  Total  Force  Level  Trajectory  For  Case  Three. 
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X  FORCE  LEVEL 


X  FORCE  LEVEL  TRAJECTORY 


o 


Figure  E.8  X  Force  Level  Trajectory  Over  Time  For  Case  Three. 
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Figure  E.9  Y  Force  Level  Trajectory  Over  lime  For  Case  Three. 
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TABLE  9 

INPUT  DATA  SET  FOR  C  S  MODEL  CASE  4 


01  05  87  TEST  OF  (1X1),  ROOTS  AT  200,200(U)  &  366, 100(S)S 

1  600  0  0  1.0  1.0  1.0  0.040  0.000  0.000  28  0  0 

0.06000  0.00000  0.00000 
0.OO000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.O00O0  0.00000  0.00000 

1  600  0  0  1.0  1.0  1.0  0.010  0.000  0.000  12  0  0 

0.01000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

345215789. DO  45635761. DO  89342761. DO 


TABLE  10 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  4 


01  05  87  TEST 
99 

1  600.  0. 
0.04972  0.00000 
0.00000  0.00000 
0.00000  0.00000 
0.00017  0.00000 
O.OUOOO  0.00000 
0.00000  0.00000 
1  600.  0. 
0.00796  0.00000 
0.00000  0.00000 
0.00000  0.00000 
0.00017  0.0O000 
0.00000  0.00000 
0.00000  0.00000 


OF  (1X1),  ROOTS  AT  200,200(U)  &  366,100(S)S 

0.  0.040  0.000  0.000  28  0  0 

0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 

0.  0.010  0.000  0.000  12  0  0 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 


Figure  E.10  Total  Force  Level  Trajectory  For  Case  Four. 


X  FORCE  LEVEL  TRAJECTORY 
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TABLE  12 

INPUT  DATA  SET  FOR  M  W  MODEL  CASE  5 


02  02  87  TEST  OF  (1X2),  CASE  5 
90 

1  500.  000.  0.  0.030  0.000  0.000  38  0  0 

0.12243  0.04733  0.00000 

0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00038  0.00032  0.00000 
0.00000  0.00000  0.00000 
0.00000  0.00000  0.00000 

2  250.300.  0.  0.050  0.100  0.000  44  31  0 

0.06156  0.00000  0.00000 

0. 04530  0.00000  0.00000 
0.00000  0.00000  0.00000 
0.00015  0.00000  0.00000 
0.00023  0.00000  0.00000 
0.00000  0.00000  0.00000 


Figure  E.  13  Total  Force  Lee  el  Trajcu.ve  i 
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Figure  E.16  XI  Force  Level  Trajectory  Over  Time  For  Case  Five. 


Figure  E.17  Y1  Force  Level  Trajectory  Over  Time  For  Case  Five. 
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Figure  E.18  Y2  Force  Level  Trajectory  Over  Time  For  Case  Five. 


TABLE  13 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  6 


01/31/87  TEST  OF  (2X2),  CASE  6,  NO  ROOTS 
99 

2  200  300  0  1.0  1.0  1.0  0.010  0.020  0.000  11 

0.01000  0.02000  0.00000 
0.01000  0.01000  0.00000 
0.00000  0.00000  0.00000 
0.00008  0.00010  0.00000 
0.000 10  0.00008  0.00000 
0.00000  0.00000  0.00000 

2  300  400  0  1.0  1.0  1.0  0.005  0.010  0.000  15 

0.02000  0.02000  0.00000 
0.02000  0.01000  0.00000 
0.00000  0.00000  0.00000 
0.00010  0.00005  0.00000 
0.00005  0.00010  0.00000 
0.00000  0.00000  0.00000 

345215789. DO  45635761. DO  89342761. DO 


16 


15 


0 


0 


TABLE  14 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  6 


01  31/87  TEST  OF  (2X2),  CASE  6,  NO  ROOTS 
99 


2  200 
0.00935  0 
0.00902  0 
0.00000  0 
0.00007  0. 
0.00009  0. 
0.00000  0. 

2  300 
0.01671  0. 
0.01638  0. 
0.00000  0. 
0.00009  0. 
0.00004  0. 
0.00000  0. 


.  300.  0.  0.010 
01828  0.00000 
00925  0.00000 
00000  0.00000 
00009  0.00000 
00007  0.00000 
00000  0.00000 
.  400.  0.  0.005 
01683  0.00000 
00823  0.00000 
00000  0.00000 
00004  0.00000 
00008  0.00000 
00000  0.00000 


0.020  0.000  11  16  0 


0.010  0.000  15  15  0 


TOTAL  Y  FORCE  LEVEL 

360  600  650  700  750  800  650  900  950  1000 


Figure  E.21  Y  Force  Level  Trajectory  Over  Time  For  Case  Six. 


XI  FORCE  LEVEL 

100  ISO  200  2S0  300  350  400  450  500 


X2  FORCE  LEVEL 
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Y1  FORCE  LEVEL  TRAJECTORY 


0  S  10  15  20  25  30  35  40  45  50  55  60  65  70  75  60  65  90  96100 


TIME 

CASE  6  (2X2)  NO  HOOTS 


Figure  E.24  Y1  Force  Level  Trajectory  Over  Time  For  Case  Six. 
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Y2  FORCE  LEVEL  TRAJECTORY 


TABLE  15 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  7 


01/31/87  TEST  OF  (2X2),  CASE  7,  2  ROOTS 
99 

2  200  300  0  1.0  1.0  1.0  0.010  0.020  0.000  11  16  0 
0.01000  0.02000  0.00000 
0.01000  0.01000  0.00000 
0.00000  0.00000  0.00000 
0.00008  0.00010  0.00000 
0.00010  0.00008  0.00000 
0.00000  0.00000  0.00000 

2  300  400  0  1.0  1.0  1.0  0.005  0.010  0.000  9  12  0 
0.02000  0.02000  0.00000 
0.02000  0.01000  0.00000 
0.00000  0.00000  0.00000 
0.00010  0.00005  0.00000 
0.00005  0.00010  0.00000 
0.00000  0.00000  0.00000 

345215789. DO  45635761.D0  89342761. DO 


TABLE  16 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  7 


01/31/87  TEST 
99 

2  200.  300. 
0.00944  0.01891 
0.00942  0.00925 
0.00000  0.00000 
0.00007  0.00009 
0.00009  0.00007 
0.00000  0.00000 
2  300.  400. 
0.01900  0.01905 
0.01821  0.00957 
0.00000  0.00000 
0.00009  0.00005 
0.00005  0.00009 
0.00000  0.00000 


OF  (2X2),  CASE  7,  2  ROOTS 

0.  0.010  0.020  0.000  11  16 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 

0.  0.005  0.010  0.000  9  12 

0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 


0 


0 


TOTAL  Y  FORCE  LEVEL 

200  290  300  390  400  430  900  930  800  830  700  750  800 


TOTAL  FORCE  LEVEL  TRAJECTORY 


CASE  7  (2X2)  2  ROOTS 


Figure  E.26  Total  Force  Level  Trajectory  For  Case  Seven. 
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X  FORCE  LEVEL  TRAJECTORY 


0  ft  10  IS  <0  25  30  35  40  45  50  55  60  65  70  75  60  65  90  95100 


TIME 

CASE  7  (2X2)  2  ROOTS 


Y  FORCE  LEVEL  TRAJECTORY 


o 


TIME 


CASE  7  (2X2)  2  ROOTS 


Figure  E.28  Y  Force  Level  Trajectory  Over  Time  For  Case  Seven. 


X2  FORCE  LEVEL  TRAJECTORY 


0  5  10  15  20  25  90  35  40  45  50  55  50  55  70  75  50  65  00  05100 

TIME 

CASE  7  (2X2)  2  ROOTS 


Figure  E.30  X2  Force  Level  Trajectory  Over  Time  For  Case  Seven. 
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100 


0  5  10  IS  to  to  SO  35  40  45  50  55  00  05  70  75  00  05  00  05 100 

TIME 

CASE  7  (2X2)  2  ROOTS 


Figure  E.31  Yl  Force  Level  Trajectory  Over  Time  For  Case  Seven. 


Figure  E.32  Y2  Force  Level  Trajectory  Over  Time  For  Case  Seven. 
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TABLE  17 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  8 


02/01/87  TEST  OF  (2X2),  CASE  8;  1  STABLE  ROOT 
99 

2  200  200  0  1.0  1.0  1.0  0.100  0.120  0.000  51  85  0 

0.04000  0.06000  0.00000 
0.05000  0.07000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00030  0.00000 
0.00010  0.00040  0.00000 
0.00000  0.00000  0.00000 

2  300  300  0  1.0  1.0  1.0  0.100  0.150  0.000  71  76  0 

0.06000  0.04000  0.00000 
0.05000  0.03000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00020  0.00000 
0.00030  0.00020  0.00000 
0.00000  0.00000  0.00000 

345215789. DO  45635761.D0  89342761. DO 


TABLE  18 

INPUT  DATA  SET  FOR  M,W  MODEL  CASE  8 


02/01/87  TEST  OF  (2X2),  CASE  8;  1  STABLE  ROOT 
99 


2  200 
0.03003  0. 
0.03652  0. 
0.00000  0. 
0.00015  0. 
0.00007  0. 
0.00000  0. 

2  300 
0.04653  0 
0.03147  0. 
0.00000  0 
0.00015  0 
0.00019  0 
0.00000  0 


.  200.  0.  0.100 
04312  0.00000 
04845  0.00000 
00000  0.00000 
00021  0.00000 
00028  0.00000 
00000  0.00000 
.  300.  0.  0.100 
03115  0.00000 
01856  0.00000 
00000  0.00000 
00016  0.00000 
00013  0.00000 
00000  0.00000 


0.120  0.000  51  85  0 


0.150  0.000  71  76  0 


TOTAL  Y  FORCE  LEVEL 

650  600  660  700  760  800  850  000  050  1000 


TOTAL  FORCE  LEVEL  TRAJECTORY 


CASE  8  (2X2)  1  STABLE  ROOT 


Figure  E.33  Total  Force  Level  Trajectory  For  Case  Eight. 
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Y  FORCE  LEVEL  TRAJECTORY 


O  5  10  15  20  25  30  35  40  45  50  55  60  65  70  75  80  05  90  95100 


TIME 

CASE  8  (2X2)  1  STABLE  ROOT 

Figure  E.35  Y  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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XI  FORCE  LEVEL  TRAJECTORY 


O 

o 

trt 


X2  FORCE  LEVEL  TRAJECTORY 


TIME 


CASE  0  (2X2)  1  STABLE  ROOT 


Figure  E.37  X2  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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Y1  FORCE  LEVEL  TRAJECTORY 


o 


TIME 


CASE  8  (2X2)  1  STABLE  ROOT 


Figure  E.38  Y1  Force  Level  Trajectory  Over  Time  For  Case  Eight. 
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TABLE  19 

INPUT  DATA  SET  FOR  C/S  MODEL  CASE  9 


02,01/87  TEST  OF  (2X2),  CASE  9,  1  UNSTABLE 
99 

2  400  400  0  1.0  1.0  1.0  0.050  0.030  0.000 
0.0S000  0.11000  0.00000 
0.12000  0.07000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00030  0.00000 
0.00020  0.00020  0.00000 
0.00000  0.00000  0.0<X>00 

2  400  400  0  1.0  1.0  1.0  0.070  0.020  0.000 
0.10000  0.05000  0.00000 
0.08000  0.15000  0.00000 
0.00000  0.00000  0.00000 
0.00020  0.00010  0.00000 
0.000 10  0.00030  0.00000 
0.00000  0.00000  0.00000 

345215789. DO  45635761. DO  89342761. DO 


TABLE  20 

INPUT  DATA  SET  FOR  M/W  MODEL  CASE  9 

02,01/87  TEST  OF  (2X2),  CASE  9,  1  UNSTABLE  ROOT 
99 

2  400.400.  0.  0.050  0.030  0.000  55  60  0 

0.04814  0.06804  0.00000 
0.07278  0.04330  0.00000 
0.00000  0.00000  0.00000 
0.00012  0.00019  0.00000 
0.00012  0.00012  0.00000 
0.00000  0.00000  0.00000 

2  400.  400.  0.  0.070  0.020  0.000  65  54  0 

0.06706  0.03417  0.00000 
0.05404  0.10125  0.00000 
0.00000  0.00000  0.00000 
0.00013  0.00007  0.00000 
0.00007  0.00020  0.00000 
0.00000  0.00000  0.00000 
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TOTAL  Y  FORCE  LEVEL 

150  250  350  450  550  050  750  050  050 


TOTAL  FORCE  LEVEL  TRAJECTORY 


TOTAL  X  FORCE  LEVEL 


CASE  9  (2X2)  1  UNSTABLE  ROOT 


Y  FORCE  LEVEL  TRAJECTORY 


TIME 

CASE  9  (2X2)  1  UNSTABLE  ROOT 


Figure  E.42  Y  Force  Level  Trajectory  Over  Time  For  Case  Nine. 
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Y1  FORCE  LEVEL  TRAJECTORY 
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Y2  FORCE  LEVEL  TRAJECTORY 


TABLE  21 

INPUT  DATA  SET  FOR  C  S  MODEL  CASE  10 


02.02/87  TEST  OF  (3X3) 

99 

3  200  150  175  1.0  1.0  1.0  0.005  0.010  0.020 
0.06000  0.02000  0.02000 
0.00000  0.00000  0.00000 
0.03000  0.01000  0.02000 
0.00000  0.00000  0.00000 
0.00000  0.00030  0.00020 
0.00000  0.00000  0.00010 

3  150  125  225  1.0  1.0  1.0  0.003  0.020  0.030 
0.05000  0.01000  0.05000 
0.00000  0.00000  0.00000 
0.02000  0.01000  0.03000 
0.00000  0.000<H)  0.00000 
0.00000  0.00020  0.00030 
0.00000  0.00000  0.00010 

3452157S9.D0  45635761. DO  89342761. DO 


4 


4 


4  11 


5  14 


TABLE  22 

INPUT  DATA  SET  FOR  M,\V  MODEL  CASE  10 


01  26  87 
99 

3  200 
0.05836  0. 
0.00000  0. 
0.02701  0. 
0.00000  0. 
0.00000  0. 
0.00000  0. 

3  150 
0.04884  0. 
O.OOXHJ  0. 
0.01884  0. 
0.00000  0 
0.00000  0 
0.00000  0 


TEST  OF  (3X3) 

.  150.  175.  0.005 
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Figure  E.47  Total  Force  Level  Trajectory-  For  Case  Ten. 
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